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0.1  Introduction:  Novel  Designs  for  Silicon-Based  Unipolar  Intersubband  Lasers 

This  contract  has  focused  on  designs  of  silicon-based  quantum  well  intersubband  lasers.  There 
are  several  strong  motivations  for  this  effort.  First,  present  infrared  lasers  are  fabricated  from 
III-V  semiconductors  and  so  require  hybrid  integration  with  silicon  based  microelectronics.  To 
have  a  folly  silicon  based  laser  monolithically  integrable  with  advanced  silicon  microelectronics 
would  be  a  significant  advance  with  respect  to  cost  and  simplicity  of  processing.  Silicon  and 
silicon/germanium  have  not  been  used  for  optical  applications  primarily  because  of  the  fact 
that  their  energy  band  gaps  are  indirect.  This  means  that  in  order  for  an  electron  to  absorb 
radiation  (a  photon)  and  be  promoted  from  the  valence  to  the  conduction  band,  crystal 
momentum  must  be  supplied  by  the  lattice  vibrations  (phonons).  Similarly,  in  order  for  an 
electron  and  hole  to  recombine  and  emit  radiation,  crystal  momentum  must  be  provided  by  the 
phonons,  substantially  reducing  the  efficiency  and  making  the  process  strongly  temperature 
dependent.  This  limitation  (the  indirect  bandgap)  can  be  avoided  by  considering  transitions 
within  the  conduction  band  or  valence  band  alone,  but  not  between  them,  i.e.  (unipolar) 
intersubband  processes  but  not  (bipolar)  interband  processes.  We  have  studied  such 
intersubband  transition  in  confined  quantum  well  structures  where  the  energy  levels  and  hence 
the  photon  energy  are  tunable  by  changing  the  dimensions  of  the  quantum  well  and  barrier. 
This  is  another  major  advantage  of  quantum  well  bandgap  engineering.  We  have  considered 
only  unipolar  hole  injection  in  the  valence  band  because  the  conduction  band  offsets  are 
negligibly  small  for  growth  on  a  pure  silicon  substrate.  A  finite  conduction  band  offset  is 
possible  for  growth  on  a  silicon/germanium  alloy  buffer  layer,  but  this  was  not  considered  in 
the  present  work  for  operation  in  the  conduction  band.  The  silicon/germanium  alloy  buffer 
layer  was  considered,  however,  for  the  strain-symmetrized  silicon/germanium  short  period 
superlattice  for  operation  in  the  valence  band,  described  later  in  this  report. 

1.1  Comparison  with  Prior  Work  and  Motivation  for  Present  Work 

There  has  been  considerable  published  work  [l]on  intersubband  lasing  using  the  III-V 
semiconductors  InGaAs  and  InAlAs.  An  injected  electron  (in  this  case)  successively  transits 
some  25  active  periods  in  an  applied  electric  field  and  has  the  possiblity  of  emitting  a  laser 
photon  in  each  period.  This  is  called  the  Quantum  Cascade  Intersubband  Laser  (QCL)  shown 
in  figure  1.  Initially  ,  when  we  considered  this  design  for  our  silicon-based  intersubband  laser 
operating  at  short  wavelengths  (1.6  microns),  it  was  apparent  that  in  order  to  energetically 
align  the  lower  lasing  level  of  a  given  period  with  the  upper  lasing  level  of  the  next  period,  an 
excessively  large  electric  field  would  be  required  ,  possibly  resulting  in  electrical  breakdown. 
To  avoid  this,  we  considered  the  Quantum  Parallel  Laser  (QPL)  in  which  the  active  periods  are 
addressed  in  parallel  rather  as  a  periodic  superlattice.  This  is  shown  in  figure  2.  Population 
inversion  has  to  be  achieved  in  the  QPL  by  bandgap  engineering  of  a  multiple  well  superlattice 
period  or  by  using  local-in-k  space  population  inversion  which  arises  from  the  nonparabolicity 
of  the  valence  subbands  (minibands  in  the  case  of  a  superlattice).  Both  cases  will  be  described. 
It  will  be  shown  that  when  coupled  wells  are  used  to  engineer  the  population  inversion,  the 
wavefonction  overlap  of  one  superlattice  period  to  the  next  is  too  small  to  give  sufficient 
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miniband  transport  for  the  QPL  to  operate  satisfactorily.  Hence,  vve  focused  on  the  local-in-k 
space  population  inversion  mechanism  which  requires  knowing  (calculating)  the  in-plane 
energy  dispersion  to  see  if  the  nonparabolicity  is  sufficient  for  this  mechanism  to  work. 


1.2  Early  models  of  Inplane  Dispersion 

Previous  work  showed  that  the  quantum  parallel  interminband  laser  in  which  each  active 
region  consists  of  coupled  quantum  wells  will  not  perform  satisfactorily.  This  is  shown  in 
figure  3.  Here  the  well  and  barrier  widths  are  adjusted  so  that  the  lifetime  of  the  upper  lasing 
level  is  longer  than  the  lifetime  of  the  lower  lasing  level.  This  can  be  engineered  because  the 
lifetimes  depend  on  the  wavefunction  overlap  which  can  be  varied  by  changing  the  width  of  the 
barrier  between  the  active  wells.  The  reason  that  this  structure  cannot  act  a  superlattice  is  that 
the  internal  barriers  within  each  active  region,  which  are  required  to  bandgap  engineer  the 
subband  lifetimes,  also  severely  decrease  the  wavefunction  overlap  from  one  active  period  to 
the  next;  this  kind  of  'lifetime  engineering'  is  required  in  silicon-based  structures  since  the 
transition  rates  are  independent  of  momentum  transfer  as  in  the  III-V's.  A  publication  resulting 
from  this  study  is  reference  2,  Photonics  Technology  Letters,  May  1997,  and  is  reproduced  in 
Appendix  A. 

A  more  tightly  coupled  superlattice  with  larger  wavefunction  overlap  and  wider  minibands  is 
needed.  One  potential  candidate  is  the  ATT  InGaAs/InAJAs  quantum  cascade  laser 
(QCL)which  does  not  achieve  'global'  population  inversion,  but  only  'local-in-k-space' 
population  inversion  [3],  A  necessary  (though  not  necessarily  sufficient)  condition  is  that  there 
be  enough  band  nonparabolicity  that  the  emitted  photon  cannot  be  reabsorbed  along  the  in¬ 
plane  dispersion  curve  except  at  kx=0.  Thus,  for  kx  ^0,  the  intersubband  energy  difference 
must  differ  from  the  photon  energy  by  more  that  the  linewidth  of  the  upper  level, see  figure  4. 
This  mechanism  requires  only  a  single  quantum  well  and  barrier  per  superlattice  period  and  so 
has  the  potential  of  providing  sufficient  interperiod  electronic  overlap  and  wide  enough 
minibands  for  coherent  transport.  The  modeling  has  to  determine  whether  the  interactions 
among  the  valence  HH  and  LH  bands  gives  rise  to  sufficient  band  nonparabolicity  for  this 
mechanism  to  work,  and  the  band  mixing  as  a  function  of  in-plane  momentum  to  obtain  the 
optical  phonon  emission  rate. 

Our  progam  was  to  study  that  mechanism  initially  for  the  GexSii.x  /Si  QPL  to  determine 
whether  there  is  sufficient  in-plane  nonparabolicity  for  this  mechanism  to  work,  and  ,  if  so,  to 
determine  the  nonparabolic  energy  dispersion  and  the  admixture  of  light  hole  and  heavy  hole 
wavefunctions  as  a  function  of  kx.  The  last  item  is  important  because  the  initial  excited  state  is 
at  kx=0  and  is  therefore  entirely  heavy-hole  like.  By  emitting  an  optical  phonon  to  a  final  state 
with  kx  *0,  the  latter  is  an  admixture  of  heavy-hole  and  light-hole  states.  Since  the  electron- 
phonon  interaction  between  states  with  different  Bloch  functions  is  small,  the  lifetime  of  the 
upper  level  will  be  increased  if  there  is  a  large  admixture  of  light-hole  wavefuction  in  the  final 
state,  enhancing  the  possibility  of  population  inversion. 
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For  the  GexSij.x  /Si  QPL,  our  procedure  was  as  follows:  To  model  the  in-plane  dispersion,  we 
used  the  6  band  (k.p)  secular  determinant  of  People  and  Sputz  [4],  This  includes  all 
interactions  in  the  valence  band  and  only  neglects  possible  interactions  with  the  conduction 
band.  Specifically,  the  secular  equation  for  the  energy  eigenvalues  is 

He[-SeLe  +  (D/2 1/2  -21/2e)2  ]  +  |a|2  [(3/2)HE  +(1/2)LE  +  SE-(D-2e)] 

+|b|2  [2Le  +Se  +2D  -4e]  +33/2  Re[  (a)2b]  =0,  (1) 

where 

He  =H-E,  Le=L-E,  Se=S-E, 

H  =Hu,  +  e,  L=  Hih  -  e,  S=1/2(L+H)-D, 

Hi!h=  -  h2/2mo  [(kx2+  ky2)(gi  +g2)+kz2(g1-2g2)], 


Hlh=  -h2/2m0  [(kx2+  ky2)(gr  g2)+kz2(gl+2g2)], 


a=3I/2h2/m0[kz(kx-iky)g3], 


b=  3 1/2h2/2m0[(kx2  -  ky2)g2-2ikxkyg3], 


D-(Hhi,- Hih),  (2) 

E  is  the  eigen-energy,  e  is  the  strain  energy,  h  is  Planck’s  constant,  the  gi  are  the  Luttinger 
inverse  effective  mass  parameters,  D  is  the  spin-orbit  splitting,  and  the  other  symbols  have  their 
usual  meanings. 


The  confinement  wavevector  kz  is  determined  from  the  solution  in  the  growth  (z)  direction 
using  appropriate  values  of  the  perpendicular  effective  masses  mz  in  the  well  and  barrier. 
Initially,  for  particular  values  of  germanium  composition  x,  zone  center  (  kx=0  )  effective 
masses  were  used  to  obtain  the  confinement  energies  and  the  corresponding  wavevectors 
kz(well)=[2mweiiE/h2]1/2  and  kz(barri er)= [2 mbnnier( V 0ffset-E)/h2] 1/2 .  For  these  latter  values,  in¬ 
plane  dispersion  was  then  calculated  as  a  function  of  kx.  While  the  dispersion  of  HH1  appeared 
plausible,  that  of  HH2  had  a  large  negative  hole  mass  which  seemed  unphysical.  One  possible 
error  is  the  assumption  that  mz  of  the  well  and  barrier  remain  constant  with  increasing  kx.  The 
procedure  was  then  modified  as  follows: 
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(1)  Values  of  kz(well)  and  k2(barrier)  corresponding  to  kx=0  were  taken, 

(2)  A  value  of  kx  *0  is  chosen, 

(3)  Values  of  m2(well)  and  mz(barrier)  are  computed  from  the  People  Sputz  dispersion 
relation[4]  using  the  finite  difference  approximation  to  the  second  derivative  of  E  with  respect 
to  kz  for  the  given  value  of  kx  and  the  values  of  kz(well)  and  k2(barrier)  of  step  (1), 

(4)  Confinement  energies  are  calculated  for  these  values  of  m2(well)  and  m2(barrier)  and 
corrected  values  of  kz(well)  and  k2(barrier)  are  obtained, 

(5)  Steps  (3)  and  (4)  are  iterated  to  convergence, 

(6)  The  energy  is  then  E=E(kx,kz),  where  kx  is  the  value  chosen  in  step(2)  and  kz  is  obtained 
from  step  (5). 

Though  producing  some  changes  in  the  in-plane  dispersion,  this  procedure  did  not  correct  the 
difficulties  described  above. 

1.3  Studies  of  GaAs/AlyGai.yAs  Quantum  Wells 

In  view  of  the  above  difficulties  in  fitting  the  in-plane  dispersion  of  GexSii-x  /Si  quantum  wells, 
we  tried  to  obtain  the  in-plane  dispersion  of  GaAs/AlyGai.yAs  quantum  wells.  This  is  a  better 
studied  material  system,  and  results  by  other  investigators  are  available.  Here  we  took 
gx=6.85,  g2=2.1,  and  g3  =2.9.  The  spin-orbit  splitting  of  GaAs  is  A=340  meV.  These 
quantities  were  assumed  to  be  the  same  for  the  barrier.  Three  cases  were  simulated:  (1)  for  the 
values  of  the  gj  cited  above,  (2)  using  the  average  of  g2  and  g3;  this  makes  the  in-plane 
dispersion  isotropic  but  does  not  give  the  correct  confinement  energies  (3)  neglecting  the  small 
difference  between  g2  and  g3,  effectively  replacing  g3  by  g2 .  This  gives  isotropy  and  also  the 
correct  confinement  energies.  However,  the  same  difficulties  were  found. 

1.4  Bastard’s  Method  for  Calculating  the  Inplane  Dispersion 

An  alternate  and  more  successful  procedure  for  calculating  the  in-plane  dispersion  is  due  to 
Bastard  [5],  It  starts  with  the  4x4  Hf8  Hamiltonian  matrix.  The  heavy-hole  and  light-hole 
envelope  wave  functions  for  finite  kx  are  expanded  in  the  envelope  wavefiinctions  at  kx=0.  The 
dimension  of  this  matrix  is  [2(m+n)x2(m+n)],  where  m  is  the  number  of  bound  heavy-hole 
(HH)  states  and  n  is  the  number  of  bound  light-hole  (LH)  states  at  kx=0.  The  matrix  elements 
are  expressed  in  terms  of  the  Luttinger  parameters  (assumed  to  be  the  same  for  the  well  and 
barrier),  the  subband  energies  at  kx=0,  and  overlap  integrals  and  momentum  matrix  elements 
between  the  HH  and  LH  envelope  wavefiinctions  at  kx=0.  The  diagonalization  of  this  matrix 
gives  the  in-plane  dispersion  as  a  function  of  kx.  We  studied  the  case  of  GexSii.x  /Si.  For 
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x=0.20,  there  are  m=2  bound  heavy-hole  levels  and  only  n— 1  light-hole  level,  while  for  both 
x=0.25  and  x=0.50,  there  are  m=2  heavy-hole  levels  and  n=2  light-hole  levels. 

We  focused  on  the  case  x=0.50,  since  this  choice  leads  to  emission  at  shorter  wavelengths. 
The  latter  case  requires  finding  the  eigenvalues  of  an  8X8  matrix  which  was  easily  done  on 
MATHCAD.  The  results  appear  to  be  physically  correct.  The  coupling  pushes  the  HH1  curve 
to  lower  energies,  while  the  HH2  curve  is  pushed  to  higher  energies,  thus  increasing  the  HH2- 
HH1  energy  difference.  We  consider  this  to  be  a  more  satisfactory  situation  than  the 
InGaAs/AlInAs  local  in  k-space  intersubband  laser,  for  which  the  conduction  band 
nonparabolicity  decreased  the  intersubband  energy  difference  with  increasing  kx.  The  amount 
of  nonparabolicity  appears  to  be  sufficient  for  this  mechanism  to  be  viable. 

We  considered  a  QPL  each  period  of  which  consists  of  25  A  of  Geo.sSio.s  (well)  and  15  A  of  Si 
(barrier),  so  the  period  is  P=40  A.  The  critical  thickness  of  the  Geo.sSio.s  layers  is  45  A  for 
stable  strain  and  100  A  for  metastable  strain;  for  the  superlattice  stack  as  a  whole, 
corresponding  to  the  average  germanium  content  of  x(avg)=0.313,  the  critical  thicknesses  are 
80  A  and  400  A  for  stable  and  metastable  strain,  respectively.  So,  with  care,  a  suitable 
pseudomorphic  superlattice  of  four  or  five  periods  can  be  grown.  The  heavy-hole  (HH)  valence 
band  offset  is  450  meV  and  the  light-hole  (LH)  offset  is  371  meV;  the  effective  masses  at  k=0 
are  nimi  =0.240  m0  and  mLH  =0.178  mo  for  the  well  and  niHH  =0.291  m0  and  iulh  =0.200  m0 
for  the  barrier.  We  solved  for  the  energy  levels  of  the  isolated  wells  and  then  for  the  minbands 
of  the  superlattice.  The  lower  HH  miniband  extends  from  96.2  meV  (at  ksi  =0)to  117  meV  (at 
kd  =tc/P),  so  its  minibandwidth  is  20.8  meV.  The  minimum  of  the  upper  HH  minband  is  at 
347.5  meV  at  the  minizone  boundary  (ks!  =7t/P);  it  is  115  meV  wide  and  extends  into  the 
continuum  above  the  HH  barrier  and  has  a  maximum  at  ks)  =0. 

Thus  the  tunnel  injected  hole  relaxes  to  the  bottom  of  the  second  miniband  at  the  minizone 
boundary.  Using  the  envelope  wavefiinctions  for  HH1,HH2,  LH1  and  LH2  at  ksi  =7t/P,  we 
computed  the  overlap  and  momentum  matrix  elements  need  to  obtain  the  in-plane  dispersion. 
These  are  of  the  form 

Ohhi.lhi  -Jdz  fnHi(z)  fuu(z),  OiiH2,  LH2  -Jdz  fnH2(z)  fu-I2(z), 

Phhi,  LH2=Jdz  fkm(z)  [d/dzfui2(z)],  (3) 


P hh2,  lh  i=Jdz  fmi2(z)  [d/dzfu-n(z)], 

Here,  the  f(z)  are  the  kx=0  envelope  functions  at  ksi  =7t/P,  where  kx  is  the  in-plane  momentum. 

As  described  in  the  previous  section,  m=2  and  n=2,  so  the  dimension  of  the  matrix  determining 
the  in-plane  dispersion  is  2(m+n)x2(m+n)  =8x8.  The  matrix  appears  explicitly  in  figure  5. 
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In  Figure  5, 

CHHi.LHj  =c(kXjky)0iU!ijLUj  (i,j=l,2), 
bHHi.LHj  =b(kXjky)Pmi,LHj,  0  ,j=l,2), 

c(kx,ky)=  3 1/2  E0  [g2(kx2-ky2)-i2g3kxky],  (4) 


b(kX)ky)=2*31/2Eo(kx-iky)g3(-i)> 

E0=h2/2m0  =3810  meV*A2 

Starting  with  the  energies  at  kx=0  ,  the  in-plane  energy  dispersion  is  plotted  as  a  function  of  kx 
for  the  HH1,  HH2,  LH1  and  LH2  levels  in  figure  6.  Also  shown  are  the  dispersion  curves  in 
the  absence  of  interaction,  denoted  by  'O'.  It  is  noted  that  the  HH1  energy  is  pushed  to  lower 
hole  energy,  increasing  its  effective  mass  in  the  growth  (z)  direction,  a  general  result  as  pointed 
out  by  Bastard.  On  the  other  hand,  the  HH2  curve  is  pushed  to  higher  energies  by  the 
interaction.  Thus  the  energy  difference  between  HH1  and  HH2  increased  with  kx  ,  a  more 
satisfactory  kind  of  nonparabolicity  than  in  the  III-V  InGaAs/AlGaAs  where  the  difference 
decreased  with  kx.  The  amount  of  nonparabolicity  predicted  by  the  above  calulation  appears  to 
be  sufficient  for  the  local-in-k-space  mechanism  to  work,  though  the  furthur  details  of  this 
mechanism  have  yet  to  be  worked  out.  A  publication  of  this  work  appeared  in  Journal  of 
Applied  Physics,  reference  [6] ,  and  is  reproduced  in  Appendix  B. 
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1.5  The  Method  of  Andreani  et  al 


The  Quantum  Parallel  laser  in  the  SiGe  material  system  operates  between  two  heavy  hole 
minbands  in  the  valence  band.  Due  to  the  k.p  interactions  among  the  valence  bands  (heavy- 
hole,  hh;  light-hole,  lh;  split-off,  so)  the  inplane  energy  dispersion  E(k//),  where  k//  is  the 
wavevector  in  the  plane  of  the  layer,  is  highly  nonparabolic.  Then,  because  the  photon  emitted 
at  k//=0  cannot  be  reabsorbed  except  at  k//=0  (due  to  the  above  nonparabolicity),  the  population 
inversion  is  local  in  k//  space  (at  kr=0) ,  rather  than  globally  over  the  entire  E(k//)  dispersion 
curve. 

Bastard’s  method  for  calculating  the  inplane  dispersion  described  in  section  1.4  assumes  that 
the  Luttinger  parameters  of  the  well  and  barrier  are  the  same,  and  ,  if  they  differ,  an  average  is 
taken.  Thus,  his  method  is  only  approximate.  Instead  of  the  GexSii.x  /Si  superlattice  on  a  Si 
substrate,  we  next  considered  the  Si/Ge  strain-symmetrized  short  period  superlattice  on  a 
GeySii.y  buffer  layer  (y=0.5  for  equal  thicknesses  of  Si  and  Ge).  This  has  a  number  of 
advantages  over  the  previously  described  Geo.sSio.s  /Si  superlattice  on  a  Si  substrate:  First,  the 
strains  in  the  Si  andGe  layers  are  opposite  in  sign  so  that  no  integrated  strain  energy  is 
accumulated  for  a  multilayered  system  and  ideally  there  is  no  limit  on  the  height  of  the 
superlattice,  as  long  as  the  individual  layer  thicknesses  are  below  their  critical  thichness  for 
pseudomorphic  growth;  this  in  not  the  case  for  the 

Geo.sSio.s  /Si  superlattice  on  a  Si  substrate  where  a  net  strain  energy  is  accumulated. 
Secondly,  compositional  rounding  of  the  Ge  concentration  in  the  growth  of  the  Geo.sSio.5  wells 
is  avoided.  This  rounding  distorts  the  shape  of  the  quantum  wellls,  changing  the  subband 
energies.  Thirdly,  alloy  scattering  is  avoided,  leading  to  longer  lifetimes.  So,  we  considered  the 
Si/Ge  strain-symmetrized  short  period  superlattice  on  a  GeySii.y  buffer  layer  .  Here  it  was 
recognized  that  the  Luttinger  parameters  (and  other  properties)  of  Si  and  Ge  are  very  different 
,  so  another  method  is  required  which  takes  these  differences  explicitly  into  account.  Such  a 
method  will  be  especially  important  for  future  work  on  e.g.  the  Si/ZnS  system  where  the  two 
materials  have  very  different  properties.  Such  a  procedure  is  given  by  Andreani,  Pasquarello, 
and  Bassani  [  7  ]  .  These  authors  exactly  solve  the  multiband  effective  mass  equation  for  the 
valence  band  envelope  functions,  analogous  to  the  solution  for  a  single  parabolic  band  with 
different  effective  masses  for  the  well  and  barrier.  Imposing  continuity  of  the  multiband 
envelope  functions  and  currents  at  the  heterointerfaces  leads  to  a  16x16  determinant  which  can 
be  reduced  to  an  8x8  determinant  by  symmetry  arguments.  From  the  zeros  of  this  determinant 
the  inplane  dispersion  E(k//)  can  be  determined.  Andreani’ s  algorithm  was  programed  in 
MATHCAD  ,  in  particular  setting  up  and  evaluating  the  8x8=64  matrix  elements  in  their 
secular  determinant  for  an  isolated  quantum  well.  The  authors  pointed  out  that  their  method 
can  be  extended  from  an  isolated  quantum  well  to  a  periodic  superlattice,  but  that  here  one 
must  deal  with  the  full  16x16  determinant.  This  extension  was  done  at  a  later  stage. 
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1.6  Calculation  of  Band  Offsets  for  the  Strain-Symmetrized  Si/Ge  Superlattice 

The  conduction  and  valence  band  offsets  of  the  Ge/Si  superlattice  on  a  Si0.Ge0.5  buffer 
layer  have  been  calculated  following  the  procedure  of  Kahan,  Chi,  and  Friedman  [8],  The 
offsets  of  the  Ge  and  Si  layers  are  found  with  respect  to  the  Si0.Ge0.5  buffer  layer ,  which 
then  determines  the  offsets  of  Si  and  Ge  with  respect  to  one  another.  The  detailed  procedure 
is  as  follows: 

Three  band  model  fheaw  hole,  light  hole,  split-off  spin  orbit  band): 

For  both  the  Si  and  Ge  layers,  the  average  valence  band  energy  including  spin-orbit 
interaction  and  hydrostatic  strain  is  given  by 
EvV  =  Ev,av  +A  /3  +  av  (2-X)e, 

where  Ev,av  is  the  average  valence  band  energy  in  the  absence  of  spin-orbit  and  hydrostatic 
strain  (  Ev,av=-7.03  eV  for  Si  and  -6.35  eV  for  Ge),  A  is  the  spin-orbit  energy 
(44  and  290  meV  for  Si  and  Ge,  respectively) ,  av  is  the  hydrostaic  deformation  potential 
(2.4  6  and  1.24  eV),  and  \=2Cnl  C22,  where  C)2  and  C22  are  the  elastic  constants 
(  C,2=  0.650x10"  dyn/cm2,  0.494x10"  dyn/cm2;  C22=  0.801x10"  dyn/cm2, 

0.684x10"  dyn/cm2  for  Si  and  Ge,  respectively),  e  is  the  inplane  strain.with  respect  to  the 
Sio.sGeo.s  buffer  layer,  e=+0.02  for  Si  (tensile  strain) ,  and  -0.02  for  Ge 
(compressive  strain).  We  calculate  Ev>av  =-6952.7  meV  for  Si  and  -6284.3  meV  for  Ge. 

The  next  step  is  to  add  to  these  quantities  the  splittings  due  to  uniaxial  shear.  These  are[4] 

dE(hh)  =  e,  (hh=heavy-hole) 

dE(lh)  =  -l/2(e+A)+l/2[  e2  +  A2  -2eA]1/2  (lh=light  hole) 
dE(so)  =  -l/2(e+A)-l/2[  e2  +  A  2  -2eA]1/2  (so=  split-off) 

We  calculate  dE(hh)  =  -77.6  mev,  dE(lh)  =  +142.27  meV,  dE(so)  =-108.66  meV  for 
Si  and  dE(hh)  =  100  mev,  dE(lh)  =-24.61  meV,  dE(so)  =-365.44  meV  for  Ge.  Finally, 
Ev(i)  =  Ev'+  dE(i),  where  i=hh,lh,so.  We  obtain 

Ev(hh)=-7030.3  meV,  Ev(lh)=-68 10.43  meV,  Ev(so)=-7061.36  meV  for  Si  and 

Ev(hh)=-6184.3  meV,  Ev(lh)=-6308.91  meV,  Ev(so)=-6649.74  meV  for  Ge.  The  valence 
band  offsets  of  Ge/  Si  are  then 

V(hh)  =  Ev(hh,  Ge)-  Ev(hh,  Si)=  836  meV,  V(lh)  =  Ev(lh,  Ge)-  Ev(lh,  Si)=  501.5  meV, 
V(so)  =  Ev(so,  Ge)-  Ev(so,  Si)=  41 1.62  meV, 
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The  band  alignment  is  type  II,  holes  are  confined  to  the  Ge  layers  and  electrons  to  the  Si  layers. 


Two  band  model  (heavy  hole,  light  hole): 


Models  for  calculating  the  inplane  dispersion  use  only  the  heavy  and  light  holes,  and  the 
split-off  spin  orbit  band  is  neglected.  So,  to  be  consistent,  the  offsets  should  be  obtained 
for  the  two  band  case.  This  is  obtained  from  the  three  band  case  by  letting  the  spin-orbit 
splitting  A-»oc,  i.e.  making  the  split-off  band  infinitely  remote.  We  then 
obtain 

dE(hh)  =  -77.6  mev  and  dE(lh)  =  +77.6  meV  for  Si,  and  dE(hh)  =  100  mev, 
dE(lh)  =  -100  meV  for  Ge  ,  so  that 
Ev(hh)=-7030.3  meV,  Ev(lh)=-6875.1  meV  for  Si  and 
Ev(hh)=-6184.3  meV,  Ev(lh)=-6384.3  meV  for  Ge. 

So  for  the  two  band  model,  the  valence  band  offsets  are  V(hh)  =  Ev(hh,  Ge)-  Ev(hh,  Si) 
=  836  meV, 

V(lh)  =  Ev(lh,  Ge)-  Ev(lh,  Si)=  490.8  meV, 

V(hh)  is  the  same  and  V(lh)  not  much  different  from  the  three  band  case. 

The  results  are  shown  in  Figure  7. 

1.7  Analysis  of  the  Strain-Symmetrized  Si/Ge  Superlattice  for  Intersubband  lasing 


As  previously  discussed  ,  the  strain- symmetrized  Si/Ge  superlattice  on  a  Si0.Ge0.5  buffer 
layer  (on  a  Si  substrate)  has  a  number  of  advantages  over  the  previously  described  Si0.Ge0.5 
/Si  superlattice  quantum  parallel  laser  diode.  First,  the  strains  in  the  Si  and  Ge  layers  are 
opposite  in  sign  so  that  no  integrated  strain  energy  is  acquired  for  a  multilayered  system  and 
ideally  there  is  no  limit  on  the  height  of  the  superlattice;  this  is  not  the  case  for  the  Si0.Ge0.5 
/Si  system.  Secondly,  compositional  rounding  of  the  Ge  concentration  in  the  growth  of  the 
SiO.Ge0.5  alloy  layers  is  avoided.  This  rounding  distorts  the  shape  of  the  quantum  wells, 
changing  the  subband  energies.  Thirdly,  alloy  scattering  is  avoided,  leading  to  longer  lifetimes. 
We  are  therefore  investigating  the  strain  symmetrized  Si/Ge  system  rather  than  the  alloy 
system. 


Summary  of  Results:  A  unipolar  p-i-p  silicon-based  intersubband  laser  consisting  of  a 
symmetrically  strained  Ge-Si  superlattice  on  a  relaxed  Si0.5  Ge0.5  buffer  layer  was  modeled 
and  analysed.  The  strain-symmetrization  removes  the  limitation  on  the  size  of  the  superlattice. 
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Analysis  of  the  in-plane  energy  dispersion  shows  that  the  population  inversion  is  local-in-k- 
space.  For  a  11ML/11  ML  superlattice  (15.4  A/15.4  A),  interminiband  lasing  between  HH2 
and  HH1  is  predicted  at  X=  2.26  mm.  From  the  envelope  functions  and  material  properties, 
the  minband  lifetimes  and  laser  gain  are  calculated.  For  a  current  density  of  10000  A/cm2,  a 
gain  of  Gl=  96  cm'1  was  calculated 

In  section  (1.4),  simulations  of  a  p-i-p  coherently  strained  Si0.5  Ge0.5/Si  superlattice  laser 
were  presented.  Calculations  were  made  of  the  local-in-k-space  population  inversion  between 
the  nonparabolic  heavy-hole  valence  subbands,  HH1  and  HH2.  Lasing  was  predicted  at  X—  5.4 
mm  with  a  lifetime  difference  of  2.4  ps.  Here  we  present  a  similar  study  for  the  symmetrically 
strained  Si/Ge  superlattice  grown  on  a  Si0.5  Ge0.5  buffer  layer.  However,  this  system  has 
several  advantages  over  the  Si0.5  Ge0.5/Si  superlattice  laser.  First,  since  no  integrated  strain 
is  accumulated  in  the  structure,  there  is  no  limitation  on  the  size  of  the  superlattice  as  long  as 
the  individual  layer  are  below  their  critical  thicknesses.  In  addition,  compositional  rounding  of 
the  Ge  concentration  is  avoided  ,  the  structures  are  simpler  to  grow,  and'  alloy  scattering  is 
eliminated.  Furthur,  the  Si/Ge  superlattice  has  large  band  offsets  so  that  lasing  wavelengths  in 
the  mid  to  near  infrared  are  possible.  Specifically,  we  propose  a  superlattice  structure  with 
equal  layer  thicknesses  of  Si  and  Ge  layers  [1 1.  monlayers  (ML)/1 1  (ML) ;  15.4  AJ  15.4  A  ] 
on  a  Si0.5  Ge0.5  buffer  layer  grown  on  a  Si  substrate.  The  Ge  layers  are  under  2% 
compressive  strain  and  the  Si  layers  under  2%  tensile  strain.  By  calculating  the  band  offsets 
with  respect  to  the  substrate,  the  band  offsets  of  the  Ge  and  Si  layers  with  respect  to  one 
another  are  obtained.  The  band  alignment  is  type  II,  with  holes  tending  to  be  in  the  Ge  layers 
and  electrons  in  the  Si  layers.  For  a  unipolar  p-i-p  structure,  hole  transport  occurs  in  the  Ge 
wells  confined  by  Si  barriers.  As  shown  above,  the  heavy-hole  valence  band  offset  at  k//=0  is 
found  to  be  Vhh  =836  meV  and  the  light  hole  offset  VU[  =501  meV,  where  k#  is  the  inplane 
crystal  momentum.  The  bandedge  energy  diagram  is  shown  in  figure  7  ,  where  the  conduction 
band  offsets  are  deduced  from  earlier  work  by  Abstreiter  et  al.  This  superlattice  gives  rise  to 
two  heavy-hole  (HH)  minibands.  The  lowest  (HH1)  is  centerred  at  252  meV  above  the  heavy- 
hole  band  edge  in  the  Ge  wells  and  is  24  meV  wide,  sufficient  for  miniband  transport  rather 
than  incoherent  hopping.  The  next  higher  heavy-hole  miniband  (HH2)  is  centerred  at  930  meV 
and  is  228  meV  wide,  extending  into  the  continuum.  Between  these  is  the  light-hole  miniband 
(LH1)  extending  from  300  to  400  meV  with  respect  to  the  heavy-hole  band  edge. 

Because  of  the  interaction  of  the  minibands,  the  inplane  energy  dispersion  is  highly 
nonparabolic.  Specifically,  the  energy  difference  between  HH1  and  HH2  increases  with 
increasing  inplane  wavevector  kx.  As  first  shown  for  the.  InGaAs/InAlAs  system,  the  photon 
emitted  at  k//=0  cannot  be  reabsorbed  at  finite'  k n  if  the  energy  difference  exceeds  the  photon 
energy  by  more  than  the  line  width  of  the  upper  level.  The  population  inversion  is  not  over  the 
entire  minbands  but  is  only  local  in  k-space  (  at  k//=0  ).  It  follows  from  that  fact  that  the 
optical  phonon  mediated  intersubband  lifetime  (HH2-HH1)  is  longer  than  the  intrasubband 
lifetime  (HH1-HH1).  In  addition,  the  effective  lifetime  of  the  lower  state  is  decreased  by  the 
probability  that  the  hole  will  tunnel  to  the  collector.  Because  of  the  large  energy  difference  in 
the  present  example,  some  12  intrasubband  cascades  are  required  to  reach  the  vicinity  of  k//=0 
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and  it  may  tunnel  at  each  cascade.  Thus  the  effective  lower  state  lifetime  is  reduced 
considerably. 


Since  the  Luttinger  parameters  of  Ge  and  Si  are  quite  different,  the  inplane  energy  dispersion 
was  calculated  by  the  method  of  Andreani  et  al  [7],  Though  their  method  can  be  generalized 
for  a  superlattice,  we  consider  an  isolated  Ge  quantum  well  surrounded  by  Si  barriers.  An 
outline  of  their  procedure  is  as  follows:  Starting  from  the  4X4  Luttinger  Hamiltonian,  the  bulk 
hole  dispersion  is  found  for  the  light  and  heavy  holes.  For  a  given  value  of  the  inplane 
momentum  k//and  energy  E,  the  corresponding  wavevectors  for  the  light  and  heavy  holes,  kz= 
ki(E)  and  kh(E),  are  found.  Next,  from  the  continuity  of  the  four  component  wavefunctions 
and  currents  at  the  well  and  barrier  interface  z-L/2,  an  8X8  determinent  is  obtained  in  which 
the  previously  obtained  values  of  ki(E)  and  kh(E)  are  substituted.  The  zeros  of  the 
determinant  D(E,  k//)  give  the  eigenvalues  E  for  the  given  value  of  k//.  In  practice,  D(E,  k//)  is 
plotted  as  a  function  of  E  to  locate  the  approximate  position  of  its  zeros  and  then  a  search 
determines  where  D(E,  k//)  =0. 

Miniband  Transport 

For  the  Ge  well,  the  Luttinger  parameters  are  giw=13.4,  g2w=4.25,  g3w=5.69,  while  for  the  Si 
barrier,  gib=4.22,  g2b=0.39,  g3b=1.44.  The  Ge  wells  are  under  2%  compressive  strain  with 
respect  to  the  buffer  layer  and  the  strain  energy  is  calculated  to  be  100  meV.  The  heavy  hole 
bandedge  lies  above  that  of  the  light  hole  bandedge  at  k//=0  by  this  amount.  Similarly,  the  Si 
barriers  are  under  2%  tensile  strain  with  a  strain  energy  of  -77  meV,  the  light  hole  edge  lying 
above  that  of  the  heavy  hole  edge  by  77  meV.  Taking  these  uniaxial  splittings  into  account, 
the  heavy-hole  valence  band  offset  at  k//=0  is  found  to  be  Vhh  =836  meV  and  the  light  hole 
offset  VLh=501  meV.  The  inplane  dispersion  of  HH1  is  repelled  toward  positive  energies  by 
an  amount  sufficient  to  render  the  population  inversion  local  in  k  space.  We  consider  a  strain 
symmetrized  superlattice  with  well  (Ge)  and  barrier  (Si)  widths  of  15.4  A  .  Simulating  the 
perpendicular  transport  for  the  heavy  holes  ( mw=0.2  mo ,  mb=0.291  mo,  Vhh—  836  meV),  we 
obtain  a  lower  heavy  hole  miniband  (HH1)  extending  from  246  meV  at  kS|=0  to  270  meV  at 
the  minizone  boundary  at  ksf=7t/P,  where  P  is  the  superlattice  period.  The  minimum  of  the 
upper  heavy  hole  miniband  (HH2)  is  at  827  meV  at  ksi=7t  /P,  just  below  the  barrier,  and 
extends  into  the  continuum.  The  optical  (lasing)  transition  occurs  at  ksi=7t/P,  the  emitted 
photon  having  an  energy  of  547  meV  corresponding  to  a  wavelength  of  X—  2.26  pm.  The 
much  stronger  nonradiative  process  is  the  emission  of  a  nonpolar  optical  phonon  to  a  state  of 
finite  inplane  momentum  on  the  lower  miniband  (HH1).  From  this  state,  the  hole  cascades 
down  to  the  minimum  at  kS|=0  by  sequentially  emitting  optical  phonons.  As  previously 
discussed,  since  the  hole  cannot  reabsorb  the  emitted  photon,  it  can  only  cascade  down  or 
tunnel  to  the  collector. 
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To  calculate  the  intersubband  transition  rate  and  therefore  the  lifetime  of  the  upper  state,  we 
need  to  know  the  envelope  wavefunctions  |i>  and  |j>  at  kS|=7r/P  from  which  we  can  calculate 
the  overlap  or  interference  factor 

Gjj(qz)  =<i  |  exp(i  qz  z)[j>, 

where  qz  is  the  phonon  wavevector  in  the  z  (growth)  direction.  Then  with 

Iij  =  J  dqz  Gij(qz), 

the  intersubband  scattering  rate  due  to  emisssion  of  an  nonpolar  optical  phonon  is  given  by 
Wij  =  [mm  Do  /  47tph(hco0p)]  [n(h©op)+l]Iij 

where  mm=0.2  m0  is  the  heavy  hole  mass  in  the  well,  D0=8.7xl08  eV/cm  is  the  optical 
deformation  potential,  p=5.32  gm/cm3  is  the  mass  density  hcoop  =  37  meV  is  the  optical 
phonon  energy  of  Ge,  and 

n(coop)=[exp(hcoop  /kT)-l]'* 

is  the  number  of  phonons  of  energy  hcoop.  Although  the  intersubband  transition  is  a  large 
wavevector  transition,  we  neglect  the  dispersion  of  the  phonon  frequency  and  take  ©op  to  be  a 
constant  (Einstein  approximation);  this  should  introduce  only  a  small  error.  We  calculate  at 
room  temperature 


Wij  =1.05xl013  A  sec'1  Iij 


For  the  above  intersubband  transition,  we  find  that  I21 — 0:0629  A *,  while  for  the  intrasubband 
transition  Iu=0.1758  A'1.  Thus  the  lifetimes  are 

t2i=  W21'1  =1.515  psec,  and  tn=  Wu'1  =0.54  psec 

Thus  population  inversion  is  established  due  to  the  fact  that  the  intrasubband  scattering  rate  is 
larger  than  the  intersubband  rate.  In  addition,  there  is  another  factor  which  leads  to  a  larger 
disparity  in  these  rates,  namely,  phonon  confinement.  Specifically,  there  is  a  large  difference  is 
the  bulk  optical  phonon  frequencies  of  the  well  and  barrier  [  vop(Si)  =15.5  THz,  vop(Ge)  =  9 
THz],  and  as  a  result,  the  optical  phonons  tend  to  be  confined.  Earlier  studies[9]of  phonon 
confinement  for  the  SiGe/Si  heterosystem  showed  that  the  intersubband  rate  was  reduced  by 
approximately  a  factor  of  3  with  respect  to  bulk  phonons  while  the  intrasubband  rate  was 
increased  by  a  factor  of  approximately  2.  Applying  these  correction  factors,  we  obtain 
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t2i'  =  (W21')*1  =4.54  psec,  and  tn'  =  (Wn  )'‘  =0.27  psec 


As  stated  above,  at  each  downward  intrasubband  cascade,  the  hole  may  tunnel  to  the 
collector.  To  obtain  the  probability  that  the  hole  reaches  the  bottom  of  the  miniband,  we  need 
the  tunnel  rate  to  the  collector  which  we  consider  next. 


Tunneling  Rate 


In  this  section,  the  tunnel  rate  form  the  lower  miniband  to  the  collector  is  calculated.  The 
centroid  of  the  lower  miniband  is  at  268  meV  above  the  heavy  hole  band  edge.  The  heavy  hole 
offset  is  V=836  meV,  the  heavy  hole  mass  of  the  Ge  well  is  mw=0.2mo,  the  mass  of  the  Si 
barrier  is 

mb=  0.291mo,  and  we  choose  the  width  of  the  barrier  to  be  5  monolayers  (7  A).  Then  with 
E0=  3810  meV  A2 

k(E)=[(E/E0)  mw]1/2 ,  k(E,V)=[V-E)mb/E0]1/2 
e(E,V)=[  k(E, V)  mw/  k(E)mb  -  k(E)mb/  k(E,V)  mw], 

The  transmission  probability  is 

T(E,a,V)=|exp(-i2k(E)a)/cosh(2k(E,V)a)+i(e(E,V)/2)sinh(2k(E,V)a)|2 
For  the  above  specified  values  of  the  parameters,  we  calculate 
T(E,a,V)=0.189 

The  mean  group  velocity  of  the  lower  miniband  is  v=AE  P/2h, 

where  AE=24  meV  is  the  width  of  the  lower  miniband  and  P=30.8  A  is  the  period.  This  gives 
v=5.6xl06cm/sec. 

Choosing  10  periods,  the  time  for  a  hole  to  transit  the  structure  is 
tlr=0.55  psec 
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The  tunelling  rate  is 
Wtun  =(1/  ttr)  T(E,a,V) 
We  calculate 


Wtun  =0.35  psec'1 ,  tum  =(  Wtun)'1  =2.9  psec. 


At  each  intrasubband  cascade,  the  probability  that  the  hole  will  not  have  tunneled  to  the 
collector  is 


Wintra/  (Wintra  "h  Wtun)  thin  /  (ttun  tintra)  0.91 

If  we  take  10  cascades,  the  effective  lifetime  of  the  lower  state  is 
h  =(0.91)10ttun=1.16psec, 

while  the  previously  calculated  lifetime  of  the  upper  level  is 
tu  =4.54  psec. 

Laser  Gain 
The  gain  is 


Gl=cu1(Nu-N1), 


where  the  cross  section 
Out  =(4rcao/n0)(EL/  T)  |zU||2, 

where  ao=l/137  is  the  fine  structure  constant,  no=3.5  is  a  mean  index  of  refraction,  El=547 
meV  is  the  lasing  energy,  and  T  is  the  linewidth  of  the  upper  level  at  ksj  =  tx/P;  we  choose  T 
=10  meV.  Finally, 

zU|  =  Jdzfu(z)zfi(z) 

is  the  optical  dipole  matrix  element,  where  fu(z)  and  fi(z)  are  the  envelope  functions  at  the 
minizone  boundary  kS|  =  7t/P  of  the  upper  and  lower  levels,  respectively.  We  obtain  zui  =  1  A. 
We  comment  on  the  small  magnitude  of  zui  later.  The  population  difference  (Nu  -  N|)  between 
the  upper  and  lower  states  is  found  from  a  rate  equation  analysis  of  the  effective  two  level 
system: 
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(Nu  -  N,)  =(J/e)tu  [1- 1,  ]  /tu  ], 

where  J  is  the  current  density  and  e  is  the  electronic  charge.  We  calculate 
(Nu  -  Ni)  =3.4x1 017  cm0  and 

Gl=  96  cm*1  for  a  current  density  of  J  =10000  A  cm’2;  of  course  the  gain  Gl  scales  directly 
with  J. 


1.8  Asymmetric  Strain-Symmetrized  Ge/Si  Interminiband  Laser 

Summary:  We  present  quantum-mechanical  analysis  of  a  proposed  Si-based  Ge/Si  superlattice 
p-i-p  laser  diode  grown  upon  a  relaxed  Sio.25Geo.75  buffer  layer.  Local-in-k-space  population 
inversion  between  hh2  and  hhl  minibands  at  the  SL  minizone  boundary  is  predicted  for  the  4.6 
um  lasing  wavelength.  Although  the  2.94  nm  quantum  well  width  is  three  times  the  0.98  nm  Si 
barrier  width,  the  superlattice  has  net  zero  strain.  A  large  dipole  matrix  element  |Z21|  =  5.7  A 
is  calculated,  leading  to  an  expected  gain  of  238  cm-1  at  3  kA/cm2  injection  current  density. 

In  the  previous  section  (1.7),  we  analyzed  a  unipolar  p-i-p  strain-symmetrized  Ge/Si 
interminiband  laser  on  a  relaxed  Sio.sGeo.s  buffer  layer.  Due  to  the  strong  valence  band 
nonparabolicity,  the  population  inversion  is  local  in  k-space.  With  the  Sio.sGeo.s  buffer  layer, 
the  well  and  barrier  dimensions  are  constrained  to  be  equal  (W  =  B)  to  insure  strain 
symmetrization  .  For  a  11  ML/11  ML  (15.4A/  15.4A)  superlattice,  lasing  was  predicted  at 
A.=2.2  im.  However,  a  small  dipole  matrix  element  z2i  =  1  A  was  obtained  resulting  in  a 
concomittant  small  laser  gain  of  Gl  =  96  cm*1  at  a  current  density  of  10  kA/cm2.  In  this 
section,  an  added  degree  of  freedom  is  gained  by  considering  a  variable  composition  buffer 
layer  and  unequal  well  and  barrier  widths  while  maintaining  overall  strain-symmetrization. 
Specifically,  we  consider  a  Ge-rich  Sio.25Geo.75  buffer  and  a  B/W  =7ML/21ML  (9.8A/  29.4A) 
superlattice.  The  wider  wells  (W=3B)  results  in  a  smaller  interminiband  photon  energy,  but  a 
larger  dipole  matrix  element  is  obtained,  |  z2i|  =  5.7  A;  this  is  significant  since  Gl  is 
proportional  to  |  z2i|2  .  Here  the  lasing  wavelength  is  X=  4.6  im  and  we  calculate  a  gain  of 
GL  =  238  cm*1  for  a  current  density  of  only  3  kA/cm2. 

Miniband  Stucture  and  Lasing  Transition 

The  band  diagram  of  our  proposed  type-II  superlattice  is  shown  in  figure  8.  The 
valence  band  offsets  were  obtained  by  the  procedure  of  Kalian  et  al  [8  ]  which  in  turn  follow 
the  procedure  of  Van  de  Waal  [10  ].  For  the  Sio.25Geo.75  buffer  layer,  the  Si  barrier  layers  are 
under  3%  tensile  strain  and  the  Ge  well  layers  are  under  1%  compressive  strain.  The 
previously  described  procedure  gives  the  following  valence  band  offsets  of  the  Ge  layers  with 
respect  to  the  Si:  VHH  =  819meV,  Vlh  =  486  meV. 
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The  location  at  k//=0  (  k >/  is  the  inplane  momentum)  of  the  heavy-hole  (HH)  minibands 
is  obtained  using  Vhh  =  819  meV,  niw  =  moe  =  0.206  mo ,  itib  =  nisi  =  0.291  mo ,  where  mo 
is  the  free  electron  mass.  The  lower  HH1  miniband  extends  from  99  meV  to  121  meV,  22  meV 
wide,  which  is  sufficient  for  coherent  transport.  The  upper  HH2  minband  extends  from  389 
meV  to  487  meV,  98  meV  in  width.  The  optical  (lasing)  transition  occurs  at  the  minzone 
boundary  kS|  =7t/P  (P=W+B  is  the  superlattice  period)  and  at  k//=0.  The  energy  of  the  emitted 
photon  is  268  meV,  corresponding  to  a  wavelength  of  X=  4.6  im  . 


Nonradiative  Processes 

The  most  rapid  nonradiative  transition  is  due  to  the  emission  of  a  nonpolar  optical 

phonon 

(acoustical  phonon  emission  is  slower  by  two  to  three  orders  of  magnitude).  The  optical 
phonon  energy  of  Ge  is  li0opt=  37  meV.  Now  although  the  268  meV  energy  difference  at  k//=0 
far  exceeds  hcoopt ,  the  optical  phonon  can  be  emitted  by  making  a  transition  from  HH2  to  HH1 
at  a  finite  value  of  k,/.  The  inplane  dispersion  of  the  superlattice  at  the  minzone  boundary  ksi 
=7i/P  was  obtained  by  extending  the  method  of  Andreani  et  al  [7]  to  a  superlattice.  The  result 
is  shown  in  figure  8  .  The  departure  from  nonparabolicity  of  HH1  and  HH2  is  clearly  sufficient 
for  the  local-in-k-space  mechanism  to  be  operative.  In  addition  to  the  two  heavy-hole  bands, 
the  intervening  light-hole  band  is  shown.  As  explained  earlier,  the  light-  hole  miniband  does 
not  seriously  affect  the  transition  HH2  to  HH1. 

The  procedure  for  obtaining  the  minband  lifetimes  was  given  earlier.  For  the 
intersubband  lifetime  we  calculate  t2i"  =  1.34  psec  and  for  the  intrasubband  lifetime  tn'  =  0.72 
psec,  so  that 

t2i"  >  ta"  as  required  for  population  inversion.  Taking  into  account  the  effect  of  phonon 
confinement  as  explained  ,  we  get  t2i  =  4  psec  and  tn'  =  0.36  psec.  Holes  occupying  the 
lower  miniband  exit  to  the  collector  by  tunneling  through  a  narrow  (4  ML=5.6  A)  Si  barrier. 
Tunneling  can  occur  at  each  step  as  the  hole  cascades  by  sequentially  emitting  optical  phonons. 
A  tunneling  time  3  psec  was  calculated  and  six  cascades  are  required.  To  obtain  the  effective 
tunnel  lifetime  ,  the  tunnel  time  must  be  multiplied  by  the  probability  of  the  hole  reaching  the 
bottom  of  the  lower  miniband.  Following  this  procedure,  the  lifetimes  of  levels  2  and  1  are 
found  to  be  t2  =  4  psec  and  ti  =  1.5  psec.  We  take  a  steady  current  density  of  J=  3000 
A/cm2  and  assume  a  linewidth  of  T=  10  meV  for  the  upper  level  2  at  ksi  =tc/P.  We  calculate  a 
population  difference  AN  =N2  -  Ni  =  1 .2xl017  cm'3  and  a  laser  gain  of  Gl  =  238  cm'1. 

For  the  proposed  quantum-parallel  laser ,  the  emitter  contact  is  a  pseudomorphic  p-Six  Gei-X 
alloy  layer  with  x  chosen  such  that  it's  quasi-Fermi  level  is  just  above  the  bottom  of  the  HH2 
miniband.  To  determine  the  Ge  composition  x  of  the  emitter,  the  bandedge  of  the  strained  Six 
Gei.x  alloy  was  calculated  with  respect  to  the  Sio.25Geo.75  buffer.  The  physical  paramenters  of 
the  alloy  were  found  by  interpolation.  It  was  found  that  for  x=0.55,  the  valence  bandedge  of 
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the  emitter  alloy  is  6  meV  above  the  bottom  of  the  upper  miniband  and,  if  heavily  doped  p- 
type,  is  suitable  for  hole  injection  through  a  thin  Si  tunnel  barrier.  The  tensile  strain  is  only 
0.8%  with  respect  to  the  buffer  layer,  so  the  emitter  can  be  a  few  hundreds  of  angstroms  thick 
and  therefore  bulk-like.  Similarly,  x  of  the  pseudomorphic  p-  Six  Gei.x  collector  contact  is 
chosen  so  that  the  bottom  of  the  HH1  lower  miniband  is  just  above  the  collector's  quasi-Fermi 
level.  For  x=0.932  (nearly  pure  Ge),  the  collector  valence  band  edge  is  an  optical  phonon 
energy  (37  meV)  below  the  bottom  of  the  lower  miniband,  thereby  facilitating  hole  transfer 
from  the  lower  miniband  to  the  collector.  The  collector  is  under  0.7%  compressive  strain.  A 
perspective  view  of  the  proposed  Fabry-Perot  laser  resonator  is  shown  in  figure  9. 

Conclusion:  Assuming  that  pure  Ge  and  Si  layers  are  used  for  the  active  superlattice  of  a 
short -period  quantum-parallel  laser,  we  find  that  an  imbalance  of  tensile  and  compressive  strain 
between  well  and  barrier  adds  a  useful  degree  of  freedom  for  optimization  of  the  laser  matrix 
element.  Unequal  strain  is  produced  by  making  the  buffer  composition  different  from 
Sio.5Geo.5,  and  if  the  ratio  W/B  is  made  inverse  to  the  strain  ratio,  the  resulting  SL  has  zero 
overall  strain.  For  a  Sio.25Geo.75  buffer,  we  calculate  a  5.7x  increase  in  |z2i|,  an  8x  gain 
enhancement,  and  a  2.  lx  increase  in  lasing  wavelength. 
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The  quantum-parallel  laser:  p+-i-p  diode  showing  lasing 
between  two  valence  mini-bands  of  Sio.5Geo.5/Si  superlattice. 
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quantum  wells  of  Fig.  1.  Diagram  shows  in-plane  dispersion  of 
valence  subbands  at  superlattice  mini-zone  boundary. 
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1.11  Appendix  A:  Quantum  Parallel  Laser :  A  Unipolar  Superlattice  Interminiband  Laser 
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Quantum  Parallel  Laser:  A  Unipolar 
Superlattice  Interminiband  Laser 

L.  Friedman.  R.  A.  Soref.  Senior  Member,  IEEE,  and  G.  Sun 


Abstract — Designs  are  given  for  parallel  i  rather  than  cascide) 
intersubband  lasing  in  superlattice  structures  of  group  rv  and 
IH-V  materials.  Structures  with  population  inversion  due  to 
coupled  wells  have  minibandwidths  too  narrow  for  coherent 
transport  and  the  incoherent  interperiod  rates  are  too  small  for 
population  of  a  significant  number  of  periods.  A  single-well  period 
with  local  population  inversion  gives  operation  in  the  coherent 
regime.  .An  example  is  given  for  the  HI-V  materials.  Both  SiO-Si 
and  Si-ZnS  are  discussed  for  intersubband  lasing  at  1--0  :/m. 

Index  Terms — Coherent  transport,  GalnAs-AJInAs.  optoelec¬ 
tronics.  semiconductor  lasers,  silicon,  SiGe-Si.  Si-Zna.  superlat¬ 
tices. 

.....  1.  Introduction 

The  quantum  cascade  laser  [1]  (QCL)  has  operated  at 
middle  and  far-infrared  wavelengths.  For  fiber-optic  com¬ 
munications.  it  would  be  advantageous  to  operate  this  laser  at 
1-3-^m  wavelengths.  However,  in  the  near  infrared,  the  QCL 
suffers  from  two  problems:  the  bias  voltage  V*  is  high  i\ a  ~ 
MVp,  where  ,Vf  is  the  number  of  periods  and  V?  is  the  voltage 
dropped  in  one  period)  and  the  electric  field  F  approaches 
or  exceeds  the  breakdown  strength  ( F  ~  ^p/P-  where  P 
is  the  period  width).  To  deal  with  these  problems,  we  have 
proposed  and  analyzed  a  new.  noncascaded  intersubband  laser 
structure,  the  quotum  parallel  laser  (QPL)  whose  voltages 
and  fields  are  low  at  1-3 -iim  emission  wavelengths.  Room- 
temperature  operation  is  expected  for  a  A  ~  1  QPL  and 
the  QPL  superlattice  structure  is  simpler  in  many  respects  than 
the  QCL  structure.  For  the  quantum-well  material,  we  have 
focused  upon  silicon  because  a  silicon-based  QPL  offers  hope 
of  low-cost  integration  with  silicon  electronic  circuits. 

n.  Laser  Structure 
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nil),  where  is  the  free  electron  mass. 

As  shown  in  Fig.  1.  each  period  includes  three  QW*s  and 
three  barriers,  with  well  widths  u/3,  “tt  increasing  from 
narrow  to  wide,  giving  rise  to  subband  levels  £3  >  £2  >  Ei* 
The  subscript  3  in  wz  (for  example)  means  that  the  level- 
3  wavefunction  is  localized  mainly  in  that  well.  Bamer  b 32 
is  wider  than  621  reducing  the  wavefunction  overlap  between 
3  and  2  more  than  that  between  2  and  1.  so  that  the  3-2 
transition  rate  is  slower  than  the  2-1  rate,  insuring  population 
inversion.  A  tunnel  barrier  67-  is  placed  between  one  period 
and  the  next  period.  The  complete  superlattice  interminiband 
laser  is  illustrated  in  Fig.  2(a).  where  we  have  shown  a  p- 
doped.  wider- band-gap  emitter  for  selective  injection  of  holes 
into  the  top  miniband  3.  plus  a  p-doped  collector  layer  at  the 
far  end  of  the  intrinsic  (undoped)  SL  that  gathers  the  charge 
carriers  transmitted  from  all  three  levels.  A  flat-band  condition 
exists  across  the  SL  because  the  applied  voltage  cancels  the 
built-in  voltage  due  to  the  work  function  difference  between 
the  injector  and  collector  contacts. 


The  QPL  is  a  unipolar  superlattice  laser  in  which  each 
active  period  lases  in  parallel  with  its  neighbor,  unlike  the 
QCL  wherein  a  single  charge  carrier  traverses  \f  periods,  pro¬ 
ducing  lasing  in  series.  It  is  however  necessary  to  insure  that 
population  inversion  is  achieved  in  each  period.  Specifically, 
in  our  superlattices  (SL),  we  initially  consider  three  coupled 
wells  for  each  period,  and  three  bound  minibands  designed  to 
achieve  population  inversion  between  minibands  3  and  -  that 
are  typically  heavv-hole  valence  minibands. 
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HI.  ANALYSIS  AND  CALCULATIONS 
We  have  performed  a  rate-equation  analysis  of  Figs.  1  and 
2(a).  The  results  are  summarized  here  (the  details  will  be 
published  elsewhere).  Here.  Wz2  is  the  phonon  scattering  rate 
from  level  3  to  level  2.  while  VV3I  and  VV21  are  the  3-to- 
1  and  2-to-l  rates.  Also.  is  the  earner  tunneling  rate 
between  level  3  of  one  period  and  level  3'  of  the  next  period 
(Wny  and  WIV  represent  2-to-2'  and  1-to-L  tunneling).  In 
a  riven  period,  a  fraction  fZy  of  the  level-3  hole  current  is 
transferred  to  level  3'  of  the  next  period  (those  fractions  are 
/rv  and  fxv  for  levels  2  and  1.  respectively).  Further,  /33;  - 
WzyfiWyy  -l-W'aa-WV)’  while  f2y  =  + 
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(b) 

Ftg.  1  (a)  Band-diagram  of  the  —  one-volt-biased  p~-i-p  for  n”-i-ni  unipo- 

lar  superlauice  mtermtmband  QPL  showing  emitter,  collector.  2-2  lasing,  and 
rapid  depopulation  of  2  by  strong  2-1  phonon  scattering  <k>i  <  'ti.l.'-  '&> 
Band-diagram  of  single  we 1 1/bam er  superlattice. 

The  above  rates  apply  to  two  limiting  cases:  1)  coherent 
miniband  transport,  and  2)  incoherent  phonon-assisted  tun¬ 
neling.  The  literature  contains  several  examples  of  coherent 
miniband  transport  [l]-[4]  for  strongly  coupled  SL's.  a  further 
example  being  the  graded  SL  coupling  regions  of  the  QCL 
[1].  Because  of  the  long  intersubband  lifetimes  noted  below 
(10-30  ps),  the  thermal  broadening  of  the  levels  is  smaller  than 
their  miniband widths.  However,  fluctuations  in  layer  widths 
and  surface  roughness  will  reduce  the  coherence  length,  so 
that  case  2)  above  is  the  more  likely  one.  especially  since  the 
calculated  minibandwidths  for  our  proposed  SiGe-Si  structure 
are  small  «5  meV).  Both  cases  will  be  discussed  below. 

We  find  that  population  inversion  .\rzi  -  -V 2»  >  0  will  occur 
at  each  period  t  (where  i  =  1.  2.  •  •  • .  -V),  provided  that  two 
conditions  are  met:  W%i  —  W2y  >  ^32  (level  2  is  depleted  at  a 
rate  faster  than  it  is  populated  from  level  3),  and  fey  ^  hr 
(tunneling  2-2'  is  much  less  probable  than  tunneling  3-3'). 
When  the  _Vf  identical  periods  are  joined  together  in  a  SL, 
levels  1.  2,  3  become  minibands  l.  2,  3.  Then  the  tunnel 
rates  \Vkk>  in  the  coherent  regime  [case  1]  are  obtained  as 
Wkk'  =  2Bk/h.  where  h  is  Planck’s  constant  and  B ^  is  the 
width  in  eV  of  miniband  k .  where  k  =  1.  2.  3.  In  turn.  B t  is 
dependent  upon  the  tunnel  barrier  width  by-  For  case  (2).  the 
Wfcfc'  are  the  rates  for  acoustical  and  optical  phonon  assisted 
transitions,  proportional  to  the  wavefuncrion  overlap  between 
adjacent  periods.  Because  W22'  <  Whs**  space  charge 

buildup  will  occur  toward  the  collector,  giving  a  redistribution 
of  the  voltage  and.  if  sufficiently  large,  make  case  2)  the 
transport  mechanism,  even  for  an  ideal  structure.  The  laser 
gain  at  stage  i  is  proportional  to  1.V3,  —  .V2i).  It  is  important 
to  note  that  the  gain  is  spatially  nonuniform  in  the  Fig.  2  QPL 


because  the  inversion  decreases  progressively  from  stage  to 
stage  and  because  of  the  space  charge. 

Our  long-term  goal  is  imerminiband  lasing  at  A  =  1.55  am. 
corresponding  to  an  800-meV  intersubband  energy  difference. 
Such  an  energy  difference  in  Si  QW’s  requires  large  valence- 
band  offsets  between  Si  and  the  barrier  material  i  l  y/i  L5  eV 
in  Fig.  1):  for  example,  offsets  provided  by  cubic  ZnS  barriers 
[5],  However,  the  relevant  carrier-phonon  scattering  rates  in 
the  Si-ZnS  heterosystem,  a  mixed  polar/nonpolar  system  are 
not  known:  these  are  under  current  investigation.  Therefore, 
to  illustrate  the  operation  of  the  QPL.  we  consider  instead 
the  Si0.:3Ge0.25-Si  SL  of  [6].  We  have  used  the  procedures 
of  [6]  to  calculate  the  rates  W32.  W31  and  M';j.  Cryogenic- 
temperature  operation  of  this  SiGe-Si  QPL  is  assumed  to  avoid 
absorptive  transitions  and  thermal  population  of  minibands  2 
and  3.  The  choice  of  bT  will  determine  how  rapidly  the  gam 
falls  off  from  stage  to  stage,  via  the  tunneling  rate  \Vzy . 

For  the  SiGe  well  widths.  Si  barrier  widths,  and  effective 
masses  given  in  the  Fig.  I  caption,  we  find  that  the  minima  of 
the  heavy-hole  minibands  HH1.  HH2.  and  HH3  are  at  52.4. 
83.1.  and  151.9  meV.  Also,  we  calculate  W*32  =  0.015  ps”1. 
tt'3l  =  0.010  ps”1.  and  W21  =  0.100  ps”1.  These  correspond 
to  a  narrow. tunnel  barrier  bj  -  20  A  selected  to  achieve  strong 
tunneling.  For  case  1).  we  calculate  that  14  n-.  W22*. 
are  0.037.  0.166.  and  1.140  ps”1.  respectively,  and  that  W*u». 
life1  ■  and  W33.  decrease  exponentially  with  by-  We  also  find 
that  fey  -  0.62  and  fey  =  0.97.  As  desired.  Wr2 1  -  W22»  » 
W32  and  fey  <  fey :  so  these  parameters  satisfy  the  QPL 
population-inversion  conditions.  The  ratio  of  laser  gain  in 
stage-. V/  of  the  SL  to  the  gain  in  stage- 1  is  ifey)M- 

For  case  2),  the  interperiod  rates  Wkk»  would  have  to  be 
comparable  to  the  intraperiod  rates  Wz2,  WZi.  W21  in  order  to 
populate  a  significant  number  of  periods.  The  calculation  of  the 
interperiod  optical- phonon-assisted  transition  rates  follows  the 
procedure  of  Palmier  ei  al.  [7].  Even  for  a  very  narrow  bamer 
between  adjacent  periods  (5  A),  we  calculate  Wzy  j\VZ2  = 
0.15.  WWW 33-  =  0.03.  and  \Vlv/WZy  =  0.002  .  too  small 
for  proper  operation.  The  reason  is  that  the  barriers  within 
each  period  required  for  population  inversion  also  decrease 
the  wavefuncrion  overlap  between  adjacent  periods.  The  same 
narrow  minibandwidths  (1-3  meV)  are  tound  for  an  analogous 
GaAs-AlGaAs  structure  [81  for  the  same  reason.  Clearly  we 
want  a  more  rightly  coupled  superlattice  exhibiting  larger 
minibandwidths.  One  very  favorable  case  is  a  superlatrice 
consisting  of  a  single  quantum  well  and  bamer  which  exhibits 
intersubband  lasing  without  global  population  inversion  but 
due  to  local-in-fc-space  inversion  [91,  shown  in  Fig.  2(b). 
For  a  superlatrice  of  48  A  Gao.47lno.33As  wells  and  40 
A  Alo.4sIno.52As  barriers  with  two  mmibands.  the  upper 
and  lower  minibandwidths  are  44  and  10  meV.  respectively, 
sufficient  for  coherent  transport.  The  optical  phonon-limited 
lifetime  of  the  upper  state  is  ~n  =  1.4  ps  [91.  We  calculate 
fey  =  0.984  and  [fey)M  =  0.09  for  .V/  =  150.  comparable 
uTthe  total  number  of  layers  of  the  QCL  [t],  [9].  However, 
due  to  space  charge  buildup  in  the  lower  miniband.  a  smaller 
number  of  periods  may  be  required  in  order  to  maintain  the 
coherent  transport.  We  plan  to  investigate  the  feasibility  of 
local-in-k  space  population  inversion  in  Si  quantum  wells. 
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Returning  to  the  Si-ZnS  near-infrared  QPL.  the  intersub¬ 
band  energy  separations  will  be  many  times  <T  at  room 
temperature,  so  that  absorptive  transitions  and  thermal  pop¬ 
ulation  of  minibands  2  and  3  will  be  minimal  at  300  K.  Due 
to  the  higher  barriers  and  stronger  confinement,  the  widths  of 
the  active  regions  and  tunnel  bamers  will  be  typically  10—20 
A.  smaller  than  those  in  SiGe-Si. 

Finally,  we  compare  the  quantum  efficiency  of  the  QPL  with 
that  of  the  QCL.  The  internal  power  efficiency  >7  is  denned  as 
the  optical  power  output  of  the  laser  per  unit  of  electrical  input 
power.  We  note  that  the  laser  voltage  is  (£3  -  £i)/«  r'or  *e 
QPL  (about  one  volt)  compared  to  >/(£*  -  £i)/e  in  the  QCL. 
If  we  let  hu  denote  the  laser  photon  energy,  let  r  «1)  denote 
the  probability  that  a  carrier  in  level  3  makes  a  radiative  rather 
than  nonradiarive  transition  to  level  2,  and  let  fzy  denote  as 
before  the  probability  that  a  carrier  in  three  tunnels  to  level 
3'  of  the  next  period  rather  than  making  a  nonradiarive  or 
radiative  transition  to  level  2  or  level  1  of  the  same  period,  then 
it  can  be  shown  that  i?(QPL)  =  rhv{\  -  [/33#i*w 5/(^3  “  £1) 
while  tj(QCL)  =  rhvjiEz  -  E L).  So  the  efficiency  ratio  is 
l  —  (/33')*v‘t  ^us*  we  anmige  10  make  jzy)M  very  small, 
approaching  zero,  then  the  QPL  and  QCL  efficiencies  are 
comparable. 

IV.  Conclusion 

Numerical  simulations  have  been  performed  on  quantum- 
parallel  lasing ‘in  Sii_rGez/Si  at  cryogenic  temperatures,  and 
estimates  given  for  the  case  of  coherent  transport.  However, 
the  minibandwiths  are  smaller  than  the  linewidths,  so  interpe¬ 
riod  transport  is  by  phonon- assisted  tunneling.  However,  the 


calculated  interperiod  rates  are  too  small.  In  general  we  believe 
that  any  period  consisting  of  coupled  wells  is  not  satisfactory. 
However,  a  superlattice  of  48  A  Gao.-j7Ino.53As  wells  and  40 
A  Alo.48lno.52As  barriers  exhibiting  lasing  due  to  local-in-/;- 
space  population  inversion  (9)  operates  in  the  coherent  regime. 
Other  structures  with  larger  miniban dwidths  require  additional 
investigation.  The  general  importance  of  the  QPL  is  that  it  is 
a  structurally  simple,  unipolar,  low-voltage  intersubband  laser, 
tunable  by  design  over  a  wide  wavelength  range. 
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Simulations  of  a  room-temperature  p-i-p  coherently  strained  Si0JGe0,5/Si  superlatuce 
quantum-parallel  laser  diode  have  been  made.  Calculations  have  been  made  ot  the  oca^-m-fc -space 
population  inversion  between  the  nonparabolic  heavy-hole  valence  mmtnands.  HH-  and  HH1. 
Lasine  is  at  5.4  p.m  and  the  optical  dipole  matrix  element  is  3.7  A.  Analysis  of  radtauve-and-phonon 
scattering  between  the  “mixed-  bands  indicates  a  lifetime  difference  between  the  upper  and  lower 
states  of 14  ps.  At  an  injected  current  density  of  5000  Aten*,  a  laser  gain  of  134  cm'  ts  calculated. 
©  1998  American  Institute  of  Physics.  [S003 1-89791 98105507-8] 


I.  INTRODUCTION 

We  recently  described  the  design  of  a  new  superlauice 
interminiband  laser  called  the  quantum-parallel  laser  <QPLi.' 
This  low-voltage  laser  operates  under  flaiband  conditions, 
and  is  much  simpler  in  structure  than  the  well-:<nown 
quantum* cascade  laser  (QCL)."  In  addition,  the  QPL  is  ca¬ 
pable  of  operating  at  near-infrared  communications  wave¬ 
lengths.  as  well  as  middle- infrared  wavelengths,  unlike  the 
QCL  in  which  the  electric  field  exceeds  the  dielectric  break¬ 
down  strength3  when  the  laser  photon  energy  approaches  0.8 

eV.  .  u 

The  purpose  of  this  article  is  to  show,  by  calculation,  the 
physical  principles  of  operation  of  such  a  silicon-based  QPL. 
The  silicon  QPL  would  be  capable  of  integration  with  ad¬ 
vanced  silicon  microelectronics.  We  shall  consider  a 
Ge ojSios/Si  QPL  comprised  of  identical,  square  quantum 
welis.  Though  there  are  four  minibands,  only  two  are  active. 
The  lasing  wavelength  is  at  5.4  /xm.  Transitions  are  in  the 
valence  band  and  population  inversion  is  local  in  k.  space 
(kx  is  the  inplane  momentum)  due  to  the  strong  valence  band 
nonparabolicitv.  Lifetimes  are  furthur  enhanced  by  phonon 
confinement.5  A  gain  of  C(,=  l34cm  1  is  calculated  for  a 
current  density  of  J  —  5000  A/cnv*. 

II.  LASER  MODEL  SYSTEM 

The  band-edge  diagram  of  our  proposed,  electrically 
pumped,  unipolar.  Si-based  QPL  is  shown  in  Fig.  L  We 
consider  a  Ge.Si^/Si  superlauice  (SL)  with  .r  =  0.5  to  ob¬ 
tain  a  large  valence  band  offset  needed  to  increase  the  lasing 
energy.  As  shown  in  Fig.  1.  we  have  chosen  the  GexSii_t 
well  width  W  =  2.5  nm.  the  Si  barrier  width  5=  1.5  run.  and 
the  period  is  />  =  W-5  =  4.0  nm.  The  SL  is  grown  on  a 
(100)  silicon  substrate.  For  .t  =  0.5.  the  critical  thickness  of 
the  compressively  strained  GetSii-«  layers  is  4  nm  for  stable 
strain  and  10  nm  for  metastable  stain.  For  the  SL  stack  as  a 
whole,  with  a  mean  Ge  concentation  of  .r=*  0.313,  the  critical 
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thickness  is  -*0  nm  for  mctastablc  strain.  Hence  with  a  period 
of  4  nm.  up  to  ten  periods  may  be  grown  pseudomorphicallv; 
four  periods  will  be  used  in  our  calcuations.  Hole  injection  is 
from  a  larger  band-gap  p-SiGe  emitter  through  a  thin  Si 
injection  barrier:  the  larger  band  gap  of  the  emitter  insures 
selective  injection  into  the  bottom  of  the  upper  HH2  mini- 
band  as  discussed  below.  Holes  are  coherently  transported 
through  the  SL  on  the  upper  miniband.  Holes  that  fall  into 
the  lower  miniband  HH1  are  collected  via  tunneling  through 
a  thin  Si  barrier,  thereby  depleting  that  miniband. 

As  shown  in  Fig.  1,  the  heavy-hole  (HH)  valence  band 
offset  at  *  =  0  is  ^=450  meV  and  the  light-hole  (LH)  off¬ 
set  is  Vlh  =  371  meV.  For  the  well,  the  effective  masses  are 
mHH= 0.24m o  and  m^O-mmo,  while  for  the  barrier 
mHH=0.29lm0  and  mLH  =  0.200m0.6  Four  minibands  are 
formed:  The  HH1  minband  with  a  21  meV  bandwidth  ex¬ 
tends  from  96  to  117  meV  at  the  miniband  zone  boundary 
k.-^fP.  where  JL  is  the  superlattice  wave  vector  and  :  is 
the  growth  direction.  The  LH1  miniband  extends  from  176 
meV  at  Jk.  =  0  to  228  meV  at  kz=ir/P.  The  HH2  miniband  is 
inverted  in  the  E~k.  plane,  with  a  minimum  at  347  meV  at 

t t/P,  is  115  meV  wide,  and  extends  into  the  continuum. 
The  LH2  miniband  is  highest  in  energy  and  plays  no  role. 
We  will  comment  on  the  role  of  the  split-off  band  later  in  the 
article.  All  energies  are  referenced  to  the  well's  HH  band 
edge. 

III.  INTERSUBBAND  LASER  OPERATION 

Figure  2  shows  a  schematic  of  hole  energy  E  as  a  func¬ 
tion  of  kx  and  k.  where  .r  is  a  direction  in  the  layer  plane. 
Holes  are  selectively  injected  near  the  very  bottom  of  the 
HH2  miniband  at  k.^rr/P.  Thus  the  initial  state  for  the 
radiative  and  nonradiative  phonon  intersubband  transitions  is 
at  the  energy  minimum  of  HH2  at  k. —  tt! P  and  kx~  0. 

For  the  radiative  processes,  photon  emission  is  vertical 
at  kx- 0  as  shown  by  the  wavy  line.  Radiation  is  from 
E^dk^TTlP)  to  £HHl(*:"7r/P)’  with  [£HH2(*:=7r/P) 
-£hhi(*-  =  7r/P)]  =  230.5  meV.  corresponding  to  a  wave¬ 
length  of  \  = 5.36  Atm.  and  is  p  polarized  (along  the  ;  axis). 

©  1998  American  Institute  of  Physics 
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FIG.  t.  The  quantum -parallel  laser:  p —  i—p  diode  showing  lasing  be¬ 
tween  two  valence  minibands  of  Si0«Ge0«/Si  superlatuce.  The  heavy  line  *.s 
the  band  edge  for  heavy  hoies  iHH>. 

This  is  the  lasing  transition.  Weaker  photon  emission  from 
^hh ;(&:=7r/P)  to  ELHl(kz  =  rr/  P)  is  s  polarized  'along  the 
x  axis)  and  has  a  finite  dipole  matrix  element  due  only  to  the 
admixture  of  conduction  band  wave  functions  into  LH1.  The 
first  transition  is  easily  distinguished  from  the  second  by  its 
polarization  and  propagation  direction  and.  with  suitable  op¬ 
tical  feedback  (resonator  geometry),  is  made  the  lasing  tran¬ 
sition.  A  ridge- waveguide  laser  cavity  with  a  TM  fundamen¬ 
tal  mode  is  most  suitable  for  the  7 -polarized  radiation.  A 
bulk  silicon  substate  could  be  employed  for  the  resonator 
because  the  spatially  averaged  index  of  refraction  for  the 
superlattice  is  larger  than  that  of  silicon  (An ~0.2).  How¬ 
ever,  a  silicon-on- insulator  (SOI)  platform  is  preferred  be¬ 
cause  the  buned  SiO^  layer  gives  a  much  larger  index  differ¬ 
ence  (An— 2)  between  the  superlattice  and  the  substrate. 
The  SOI  provides  high  Q.  good  mode  confinement,  and  a 
high  mode  overlap  factor.  Figure  3  illustrates  the  proposed 
rib  laser  configuration  with  vertical  end  facets  in  a  Fabry- 
Perot  cavity. 

Consider  next  the  nonradiauve  processes:  At  the  bottom 
of  the  HH2  miniband  at  kz  —  7 t/P.  the  miniband  group  ve¬ 
locity  is  very  small.  At  room  temperature,  the  probability  of 
optical  phonon  emission  is  proportional  to  [1  -rt(.oiopij], 
where 

n(  ^opt)  =  [exp(  h  u)opt/kT)  -  1  ]  *  * ,  (1) 

E(kx,  kz» 

4 


FIG.  2.  Dispersion  of  lowest-three  valence  minibands  of  Si0jGeoj'5i  su- 
periaidce  as  a  funcuon  of  kt  and  k. ,  W;i  is  the  imersuband  transition  rare 
and  IV. ,  is  the  inuasuband  rate. 


FIG.  3.  Wavegtiided  unipolar  SiGe/Si  •iuperlatuce  .mermimband  laser  :n 
silicon-on-insulator.  Fabry-Peroi  resonator  :s  shown. 


The  emission  probability  is  a  factor  of  7.25  larger  than  that 
of  phonon  absorption,  which  is  proportional  to  n(  o>opi).  and 
this  factor  is  still  larger  at  lower  temperatures.  For  the 
Ge05Si#)  j  alloy,  the  averaged  optical  phonon  energy  is  taken 
as  51  meV.  Holes  absorbing  a  phonon  are  scattered  to  the 
center  of  the  HH2  mmiband  where  they  rapidly  transit  the 
superlatticc  and  are  not  available  for  lasing.  However,  at  300 
K.  the  probability  of  thermal  promotion  of  holes  to  the  center 
of  the  upper  miniband  via  optical  phonon  absorption  is  only 
0.12.  By  careful  selective  injection  into  the  bottom  of  the 
upper  miniband  state  with  small  group  velocity,  hole  leakage 
current  to  the  p-type  Si  collector  due  10  tunneling  from  this 
excited  state  is  avoided,  to  within  the  above  phonon  branch¬ 
ing  ratio.  Another  reason  for  believing  that  300  K  laser  op¬ 
eration  is  feasible  is  that  the  thermal  excitation  of  holes  from 
HHI  to  HH2  is  negligible  at  300  K. 

Turning  now  to  the  nonradiative  HH2-HH1  transition  in 
Fig.  2.  emission  of  an  optical  phonon  from  HH2  occurs  to 
superlattice  states  of  finite  in-planc  momentum,  kx.  from 
which  the  holes  cascade  down  via  intrasubband  phonon 
emission.  More  exactly,  as  shown,  the  final  states  are  on  a 
curve  which  is  the  intersection  of  the  energy  surface 

horizontal  plane 

=  =  ft  a>0Pf  We  will  assume  that  the 

holes  relaxing  into  HHI  are  distributed  over  the  lower  mini¬ 
band  and  transit  the  SL  at  the  mean  miniband  velocity  in 
calculating  the  tunnel  time  below.  However,  in  the  calcula¬ 
tion  of  optical  phonon  scattering  rates,  we  shall  consider 
only  the  in-plane  transition  at  kz~Tr/P,  in  the  interest  of 
calculational  simplicity. 

Due  to  the  interactions  between  the  HH  and  LH  states, 
the  in-plane  valence  band  dispersion  is  highly  nonparabolic. 
As  shown  in  Fig.  A  £HH1(fct)  is  forced  to  lower  hole  ener¬ 
gies.  £hh:(^t)  *s  forced  to  higher  hole  energies,  and 
£ LHif kx)  becomes  electron-likc  at  £t*=0.  Thus  the  energy 
difference  [£hhi(*t)  increases  with  increasing 

kx  and  the  photon  emitted  at  kx~  0  cannot  be  reabsorbed  as 
the  hole  cascades  down  the  £hhi(^i)  dispersion  curve,  until 
it  reaches  the  minimum.  More  exactly,  reabsorption  will  oc¬ 
cur  when  the  absorbed  photon  is  within  the  linewidth  of  the 
upper  level.  We  therefore  have  a  local  in  A* -space  inversion 
as  had  been  described  for  the  InGaAs/lnAlAs  system  by 
Faist  et  ai /  except  that  here  the  intersubband  energy  differ- 
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FIG  i.  In-plane  dispersion  E=Elkx>  of  HH1.  HH2.  and  LHl  bands  ai 
mimzone  boundarv  i.-ir/P  of  Si^Geoj/Si  suptrlamce.  sol.d  lmev 
Dashed  lines  are  in  the  absence  of  interaction  i(.  =  c  =  0).  The  .one  -ashed 
straight  line  shows  the  intereubband  optical  phonon  emission.  Decreasing 
electron  energy  (increasing  hole  energy!  is  in  the  downward  arecaon. 


ence  increases  rather  than  decreases  with  kx :  expiictt  calcu¬ 
lation  shows  that  at  fct=0.02A'\  [£hhi(*t) ’"^hkc^'WJ 
has  increased  by  50  meV  as  compared  to  its  value  at  kg 
=  0.  At  each  step  of  the  cascade,  the  hole  can  either  make  an 
downward  intrasubband  transition  by  optical  phonon  emis¬ 
sion  or  tunnel  to  the  collector.  Thus  the  tunnel  lifetime  of  the 
lower  level  is  decreased  by  the  probability  that  the  hole  cas¬ 
cades  to  the  minimum  without  tunneling  to  the  collector. 
This  decreases  the  effective  lifetime  of  the  lower  level  and  so 
increases  the  lifetime  difference  between  the  upper  and  lower 
lasing  levels,  enhancing  the  population  inversion. 

This  lifetime  difference  is  further  enhanced  by  phonon 
confinement:  namely,  earlier  calculations5  of  phonons  con¬ 
fined  to  a  Ge0.5Sio.5  layer  with  Si  barriers  gives  a  factor  of  3 
enhancement  of  the  intersubband  lifetimes  due  to  confine¬ 
ment  [cf.  Fig.  5(a)  of  Ref.  5],  while  the  intrasubband  lifetime 
is  reduced  approximately  by  a  factor  of  2  [cf.  Fig.  4(a)  of 

Ref.  5].  .  .  . 

It  may  appear  that  the  existence  of  two  nonradian ve  de¬ 
cay  paths  for  the  hole  in  the  initial  HH2  state,  namely  to 
HH1  and  to  LHl.  would  result  in  too  small  a  lifetime  for  the 
initial  state.  However,  with  increasing  in-plane  momentum 
kx,  the  HH1  envelope  function  acquires  an  admixture  of 
light-hole  wave  function,  while  the  LHl  acquires  heavy-hole 
character.  It  is  known  for  bulk  materials  that  the  HH  to  LH 
transition  rate  is  one  to  two  orders  of  magnitude  smaller  than 
the  HH  to  HH  rate7  due  to  the  smallness  of  the  overlap  of  the 
periodic  part  of  the  Bloch  functions  belonging  to  different 
energy  bands.3  Wc  assume  that  same  applies  to  confined 
valence-band  states  so  that  the  HH  to  LH  rate  may  be  ne¬ 
glected  compared  to  the  HH  to  HH  rate.  Thus  the  net  optical 
phonon  emission  rate  from  HH2  to  *‘HH1'’  is  decreased,  this 
being  compensated  for  by  an  increase  in  the  HH2  to  *‘LH1 
rate.  Thus  the  intersubband  scattering  rate  from  HH2  is  not 
increased  due  to  the  presence  of  LHl.  This  is  explicitly  ac¬ 
counted  for  in  the  calculations  below. 
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IV.  MODELING  OF  SiGe  IN-PLANE  DISPERSION 

The  SiGe  quantum  wells”  in-plane  dispersion  diagram  is 
required  for  the  laser  gain  calculation.  To  model  the  in-plane 
valence  band  dispersion,  we  follow  the  procedure  of 
Bastard.9  The  starting  point  is  the  [*1X4]  Luttinger  matrix 
multipWing  the  4-column  state  vector  whose  components  arc 
the  envelooe  functions  at  kx±0  corresponding  to  the  3/2.3/ 

2).  3/2.- i/2).  3/2. 1/2).  and  3/2.-3/2>  states,  respectively, 
and  is  given  by  Eqs.  (AH  and  (A2)  of  the  Appendix.  As 
shown  there  in  detail,  the  procedure  is  to  expand  the  enve¬ 
lope  functions  at  kx* 0  in  the  bound  state  envelope  finctions 
aufc.  =  0.  cf.  Eq.  (A3).  From  these  we  obtain  Eqs.  (A4MA7). 

This  procedure  leads  to  a  2(m  —  n)X2(m-nj  matrix, 
the  diagonalization  of  which  gives  the  E(kx)  diagram.  For 
the  particular  superlartice  strucrure  described  above,  there  are 
two  bound  HH  states  [m  =  2)  and  two  bound  LH  states  in 
=  2).  at  &.=*  tr/P.  As  explained  in  the  Appendix,  the  matrix 
is  of  dimension  2(m-n)X2(m->i).  leading  to  the  8X3 
shown  in  Table  I.  The  eigenvalues  of  this  matrix  as  a  func¬ 
tion  of  kx  gives  the  in-plane  dispersion. 

V.  RESULTS  OF  IN-PLANE  DISPERSION 
CALCULATION 

Because  the  wells  have  a  center  of  symmetry',  the  off- 
diasonal  components  of  0  le.g..  Ohhi.lhz)  ^  diagonal 
‘  components  of  P  [e.g..  Phhi.lhiJ  g*ven  Ecl*  (A9)‘  vanish' 
For  the  other  components,  by  numerical  integration,  we  ob¬ 
tain:  Ohhi.lhi  =  0.9997,  Ohhilhz®  0.9932.  Phhi.lh: 

=  "Plh:.hhi  =  0.0798.  Phhz.lhi=  “  ^lhi.hh:-0'07820- 
Solving  for  the  eigenvalues  of  the  matrix  by  standard  proce¬ 
dures.  we  obtain  the  in-plane  dispersion  curves  shown  in  Fig. 

4;  here  electron  energy  is  negative  in  the  downward  direction 
corresponding  to  positive  hole  energy.  The  dotted  curves  are 
the  parabolic  dispersions  in  the  absence  of  interaction  (b 
=  c  =  0).  It  is  noted  that  the  lowest  eigenvalue  £hhi(*x)  is 
repelled  to  lower  energies,  increasing  its  in-plane  ^effective 
mass,  a  general  result  as  pointed  out  by  Bastard.  On  the 
other  hand.  Em(kx)  is  repelled  to  higher  hole  energies  so 
that  the  difference  [£hhi(*t)  “  increases  with  A,; 

this  is  opposite  to  the  usual  conduction  band  nonparabolicity 
which  gives  a  larger  mass  with  increasing  energy,  thereby 
decreasing  the  energy  difference  with  kx .  This  situation  is 
consider  more  satisfactory  in  supressing  intersubbsnd  ab¬ 
sorption  at  finite  kx.  The  intermediate  £lhi(^x)  dispersion  is 
electron-like  at  kx  =  0. 

At  k.-irIP,  the  intersubband  optical  phonon  emission 
from  HH2  at  *x  =  0  goes  to  HH1  at  =  0.1485  A~!.  This  is 
180  meV  above  the  HH1  minimum  at  *,  =  0.  so  three  intra¬ 
subband  optical  phonon  emissions  occur  (hu)opl-5\  meV). 
brineins  the  hole  to  within  25  meV  from  the  minimum. 
Similarly,  the  intersubband  transition  from  HH2  goes  to  LHl 
at  **  =  0.115  A"1. 

Finally,  we  comment  on  the  neglect  of  the  split-off  (SO) 
band  in  the  above  analysis.  Even  at  *  =  0.  the  SO  band  will 
interact  with  the  LH  band10  and  perturb  the  band-edge  ener¬ 
gies  by  tens  of  meV;  this  is  neglected  here.  The  valence  band 
offset  of  the  SO  band  is  shallow,  only  1 10  meV  with  respect 
to  the  SO  level  of  the  Si  barrier.6  This  together  with  the  small 
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SO  effective  mass  *0.07  meV  at  r  =  0.5)  gives  a  wide  SOI 
miniband  whose  energy  at  kz  ~  tt/  P  is  526  me V .  above  die 
minimum  of  HH2  at  347  meV.  The  SO  band  will  generally 
repel  all  the  bands  to  higher  (more  positive)  energies/1  We 
are  assuming  that  the  altered  nonparabolicity  allows  the 
local- in-L -space  mechanism  to  be  operative.  The  inclusion  of 
the  SO  band  would  require  further  investigation. 

VI.  NUMERICAL  ESTIMATES  OF  TRANSITION  RATES 

Assuming  unconfined  phonons,  the  nonpolar  optical 
phonon  transition  rate  from  state  i  to  state  j  is  given  by-**'5 

Wfj  =  (m*Do/47rpft:u>opt)[n(  (Dopt)  —  1/2=  l/2]/y.  (2) 

where5  D0  =  8.7X10*  eV/cm  is  the  optical  de¬ 

formation  potential. 

p  =  3.83  gm/cm3  (3) 

is  the  mass  density  of  the  Geo^Sioj  alloy.  ^  “_5 1  meV  is 

the  optical  phonon  energy  of  the  alloy,*  n(<i>5pi) 

-  l/[exp(ftajopt)  - 1  ]  is  the  phonon  population.  li} 

=  /  dq^Gijiq.y1.  and  Gif{qz)  =  {i\txp{iq:z)\j)  is  the  hoie- 
phonon  overlap  factor.  By  numerical  integration,  we  obtain 
^hh2-hhi=0.115  A*1  (intersubband).  / hhi-hhi ”0.255  A" 1 
(inirasubband),  giving 

Twer”  1/WHH2-HHl”0*74  Ps* 

rintn"  1/WHHI-HH1=0-33  Ps- 

Applying  the  previously  mentioned  corrections  due  to 
phonon  confinement.5 

riaIcr=2.22  psec.  7intn=0,16  ps.  (5) 

Strictly,  these  correction  factors  apply  to  an  isolated  well3 
but  here  we  assume  they  apply  to  a  superlattice.  To  take  into 
account  the  intersubband  transitions  to  both  HH1  and  to 
LH1,  we  calculate  the  normalized  eigenvectors  correspond¬ 
ing  to  the  specific  final  state  energies.  As  stated  above,  the 
initial  HH2  state  is  at  meV.  k.  =  irlP,  kx- 0. 

The  final  HH1  state  at  £hhi  =  ~296  meV.  kz=irlP,  kx 
=  0.1485  A’1.  Substituting  these  values  into  the  eight  equa¬ 
tions  given  by  Eqs.  (A4)-(A7),  we  find  two  degenerate 
eigenvectors  corresponding  to  the  Kramers  degenerate  levels 
13/2,3/2)  and  13/2.— 3/2),  namely,  (cel,az./3l  .yz. 

<rx  i0-2)  =  (0.537,0,0.73 1,0,0. -0.287/. 0.0.3060  and  (0. 
-0.306<\ 0,0.2871,0.73 1.0.537,0).  For  the  final  LH1  state  at 
kz=irlP .  fcx  =  0.115  A~!,  we  obtain:  (0.446.0,-0.433.0.0. 
—  0.6 li.0.- 0.490  and  its  degenerate  counterpart.  Subtract¬ 
ing  out  the  probabilities  of  the  light-hole  components,  the  net 
effect  of  the  wavefunction  mixing  is  to  increase  by  a 
factor  of  approximately  (0.8) -1  =  1.25.  Thus,  'nIer=2.8  Ps  ^ 
compared  to  rsntcr=2.22  ps  given  in  Eq.  (5). 

VII.  TUNNELING  RATE  TO  COLLECTOR 

Holes  exit  the  superlattice  from  the  lower  miniband 
through  a  thin  tunnel  barrier  to  the  collector.  The  center  of 
the  lower  miniband  is  at  an  energy  of  107  meV.  the  barrier 
height  is  V  =  VHH=450  meV,  {mwlm0)~ 0.24.  (m^/m0) 
=  0.29.  and  we  take  the  Si  barrier  to  be  five  monolayers 


wide.  2a  =  (5  ML)  =0.7  nm.  Denning  the  tunnel  amplitude14 
A(Em.V)  as  the  amplitude  of  the  incident  wave  transmitted 
through  the  barrier,  we  have 

AiE.a.V)  =exp[t2fc(  £)a]/{cosh[2tf(£.  V)a] 

~t[ei  £.  V)/2)sinh(  2/rt  £,  V)a]}, 

where  fe(£)  =  [(£/£o)(mw/mo)]1'2.  *(£. V)  =  [( V-E)/E^) 
X(mj,/m0)] 1/2.  and 

€(E.V)  =  {KiE.V)/kiE)]{mw/mb)-[k{E)UiE'V)] 

X(mj,/mw).  (6) 

The  tunnel  probability  is  7(£.a.  V)-\A(E.& V)r. 

We  calculate:  T(E.a.  V)  =  0.291.  This  is  assumed  to  be 
independent  of  the  transverse  (in-plane)  energy.  The  mean 
velocity  of  the  lower  miniband  is  u  =  A £*  P/2h  =  6.4 
x  106  cm/s.  where  A£=21  meV  is  the  width  of  the  lower 
miniband.  Assuming  the  superlattice  to  consist  of  four  peri¬ 
ods.  its  length  is  £=16  nm  and  the  time  to  transit  the  struc¬ 
ture  is  r;r=  Liu  =  0.25  ps.  The  tunnel  rate  is  the  transit  rate 
multiplied  by  the  tunnel  probability  Wtun=(l/rtr)*T(£. 
a .  V) ,  and  the  tunnel  time  is 

Tun~  1/Wnjn=0.86  ps.  (7) 

As  pointed  out  by  the  Lucent  Technologies  work  on  the 
local  in  it -space  mechanism.4  at  each  stage  of  the  intrasub- 
band  cascade,  the  hole  may  tunnel  to  the  collector.  The  prob¬ 
ability  at  each  stage  that  it  cascades  down  rather  that  tunnels 
is  Wi;„/(Wintra-V^)  =  r^n/(rtun-rincra)  =  0.84.  There  are 
three  cascades,  so  the  hole  can  tunnel  or  undergo  intrasub- 
band  scattering  at  four  junctures.  Thus  the  probability  that 
the  hole  reaches  the  minimum  is  [i'tun/(rtun- 7inlJ.o)]4  =  0.5, 
and  the  effective  tunnel  time  is  (7am)ctf=rmn*[7tua/(7lun 
^unlra)\i  =  0A2  ps. 

In  summary,  having  applied  the  correction  factors,  the 
lifetimes  are 

rinier=~8  Ps*  (g) 

(-„)«"=  0.43  ps. 

VIII.  EMITTER  AND  COLLECTOR  CONTACTS 

The  collector  is  (bulk)  Si0.5Ge05  doped  p  type  with  bo¬ 
ron  to  2X  1019  cm*3.  This  brings  the  quasi-Fermi  level  of 
the  collector  30  meV  above  the  band  edge,  which  is  well 
below  the  bottom  of  the  lower  miniband  at  96  me  V.  We  have 
assumed  a  3D  density  of  states  effective  heavv-hole  mass15 
of  m=  1.0m0  at  300  K,  For  the  emitter  contact,  a  Geo.3Si0.7 
alloy  is  chosen.  Its  heavy-hole  band  edge  is  260  meV  below 
that  of  the  silicon  substrate  and  the  silicon  barriers  of  the 
superlattice.  When  doped  to  2X  1019  cm-3  boron,  the  quasi- 
Fermi  level  of  the  emitter  is  just  at  the  bottom  of  the  upper 
miniband.  as  required.  The  thickness  of  the  silicon  tunnel 
barrier  from  the  emitter  is  chosen  to  be  the  same  as  the 
collector  barrier,  5  ML  or  7  A. 
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IX.  LASER  GAIN 

The  gain  is  written1 2 

Gt*o-:i0V:-:V.)r.  \9\ 

where  the  lasing  cross  section  is  (rzl  =  {^a0fnQ) 

X ( / T) !Chh2-hhi  with  <2q  =  1/137,  the  fine  scructurc 
constant.  n0=3.5.  the  refractive  index.  EL  =  230.5  meV.  the 
lasing  energy.  V  is  the  linewidth  of  £kh:^i)  at  k.-rr/P. 
since  the  initial  state  is  at  k.=  m/P;  we  choose  T  =  10  meV. 
2hh2-hhi  *s  ^  dipole  matrix  clement  for  the  lasing  transi¬ 
tion.  where  Chh:-hhi  =  =  J  o^hh2^)*-^hk: 

x(cWc. 

Here,  :HH1)=<£hhi(;:)  and  HH2)  =  : )  are  die  en¬ 

velope  functions  at  &T=0  and  k.-ir/P  and  satisfy  antiperi- 
odic  boundary  conditions,  “<£hhi(-  =  ^)[*  =  !•-]• 

We  obtain  Chh:-hhl-3*7  A.  The  factor  r=  7/8  in  Eq.  (91  is 
the  phonon  branching  ratio  described  in  Sec.  II,  cf.  Eq.  1 1)  et 
seq. 

The  population  inversion  difference  is  found  from  a  rate 
equation  for  the  coupled  two  level  system: 

(N2“*/V  •)  —  (-//£)  ^interD  “1  *wn)<^/  "inter)  * 

where  J  is  the  charge  current  density  and  e  is  the  electronic 
charge.  We  calculate  (/V:“.V,)  =  7.4X  10‘°  cm"'.  and  G L 
«  134  cm-1  for  7=5000  A/cm2.  The  gain  is  proportional  to 
the  current  density,  and  since  the  emitter  and  collector  Fermi 
levels  are  at  the  bottom  of  the  upper  and  lower  minibands, 
respectively,  the  applied  voltage  required  for  the  fiatband 
condition  is  ^  =  0.23  V.  GL  is  inversely  proportional  to  the 
value  taken  for  f\ 

-X.  CONCLUSIONS 

Simulations  of  a  room-temperature  p-i-p  coherently 
strained  SiojGeo.j/Si  superlatticc  quantum-parallel  laser  di¬ 
ode  have  been  made  and  a  waveguided  SOI  resonator  struc¬ 
ture  has  been  proposed.  Calculations  have  been  made  of  the 
local-in-fc-space  population  inversion  between  the  nonpara¬ 
bolic  HH2  and  HH1  valence  minibands  at  the  5.4  psn  laser 
wavelength  where  the  dipole  matrix  element  is  3.7  A.  Analy¬ 
sis  of  radiative-and-phonon  scattering  between  the  “mixed” 
bands  indicates  a  laser  lifetime  difference  of  2.4  ps.  The 
non-radiative  intersubband  optical  phonon  scattering  is 
slower  than  the  in- plane  intrasubband  scattering.  The  inter¬ 
subband  lifetime  is  enhanced  and  the  intrasubband  lifetime 
decreased  by  phonon  confinement.  The  effective  lifetime  of 
the  lower  state  is  furthur  decreased  by  the  probability  of 


sequential  tunneling  to  the  collector.  At  a  current  density  of 
5000  A/cm',  a  laser  gain  of  134cm"1  is  calculated. 


APPENDIX:  MODELING  OF  THE  IN-PLANE 
DISPERSION 

The  SiGe  quantum  wells’  in-plane  dispersion  diagram  is 
required  for  the  laser  gain  calculation.  To  model  the  in-plane 
valence  band  dispersion,  we  follow  the  procedure  of 
Bastard.9  The  stoning  point  is  the  [4X4]  Luttinger  matrix 
multiplying  the  4-column  state  vector  whose  components  are 
the  envelope  functions  at  £x*0  corresponding  to  the  ‘3/2.3/ 
2),  3/2, -1/2).  3/2. 1/2).  and  ;3/2.-3/2)  states,  respectively. 
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Hi* 

0 

-b* 


b 

0 

#ih 

c* 


0 

-b 

c 

^hhi 


(Al) 


where 


IA2) 


H hh=  -  £ oi< K  ~  y\  -  ri  “  -  7:)]  ^ 

“ E0Hk;-kl)i  y,  -  ?:)“*;(  y,  -2y:)]  -  f. 
c28  c(kx  )  *  3 1/2£0[  yz(k;  -k;)-  i2  y  }kxky] . 

b  =  b( kt  ,ky)  =  2 *3  [nE0(kx -  iky)  y3(  -  i)dfdz 
=  b~  [k^.k^dldz, 

where  e  is  the  uniaxial  strain  energy.  E0  =  ft:/2m0 
-  3810  meV*  A2,  and  the  yf  (i'=  1,2,3)  are  the  Luttinger 
parameters.  At  finite  kx .,  the  envelope  functions  are  ex¬ 
panded  in  the  envelope  functions  at  /ct  =  0: 

where  there  are  m  bound  HH  states  and  n  bound  LH  states  at 
kz- 0.  Thus  there  are  a  total  of  2(m~n)  coefficients.  Equa¬ 
tions  (Al)  and  (A3)  are  substituted  into  the  eigenvalue  equa¬ 
tion  tfrg  ^  =  E'l'.  Considering  the  first  row.  multiplying  to 
the  left  by  (O*®)*  and  integrating,  we  obtain 

a, ( i \*cmnMjPj~'Zj-\  ^mnui  7j  = 


*  = 


(A3) 


:  L...JR 


(A4) 


TABLE  I.  8x8  matrix  determining  in-plane  dispersion. 
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where  *hh,lh )=($'?'' biQ'f). 

Following  the  same  procedure  for  [he  second,  third.  and 
fourth  rows,  we  obtain 

2 ;  a  LmcLHiKH;a;  ”  A(  ^LH*  I  “  0- 

l  =  1 . n:  >A5) 

Lm^LHiHH/*/  ~  r.(£LH»  ~  E)  ”  =  0. 

i=l . n:  fA6) 

“  ■  l.«^LHiHH/ A  *  1  ?/  ~  ®  0. 

i  =  l . m:  (AT) 

where  the  following  relations  hold: 

This  procedure  leads  to  a  2(m  —  fl)  X  2(m  —  n )  matrix,  the 
diagonalizauon  of  which  gives  the  E(kx)  diagram.  For  the 
particular  superlattice  structure  described  above,  there  are 
two  bound  HH  states  (m  =  2)  and  two  bound  LH  states  in 
=  2),  a t1cz  =  ir/P,  and  so  the  matnx  is  of  dimension  3*8 
and  is  shown  in  Table  I. 

In  this  matrix,  alf entries  are  functions  of  kx  and  hy.  and 
we  have  set  Jty*0.  Explicitly, 


C  HHi.  LH/  “~C\kx  .ky)0  HHi.LH;  {i.j-12),  (A8) 

where 

HHI.LH  1  ~  J 

^hh:.lh:=  J  ^9) 

^HHI.LH:”  j 
P HH2.LH1  =  j 
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